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SECTION  I 


INTRODUCTION  AND  SUMMARY 

This  ’  eport  presents  and  discusses  work  performed  and  results 
obtained  under  Contract  C-l  04-65  from  28  April  through  15  October  1966.  During 
that  period,  the  hypocenter  and  magnitude  progiams  were  tested  and  then  used 
to  process  January  1964  data  at  the  computer  facilities  of  the  Environmental 
Science  Services  Administration  (ESSA),  Suitland,  Maryland,  using  the  CDC 
6600  computer.  Results  of  this  processing  are  shown#in  Appendix  A. 

Other  results  obtained  in  the  reporting  period  are  given  in 
Section  II  and  in  the  appendixes.  Results  presented  in  Section  II  can  be 
summarized  as  follows: 

•  The  relationship  between  m^  and  restrained  to  a 

slope  of  0.63,  is  given  by  m^  =  0.63  M^  +  1.77  as 

compared  to  Richter's  (1958)  m,  =  0.  63  M  H-  2.  5 

b  s 

•  Magnitudes  comparable  to  mfa  but  based  on  long- 
period  P  amplitudes  average  approximately  0.  7 
units  higher  than  m^ 

•  Vertical-component  surface -wave  magnitudes 
average  about  0.  5  units  higher  than  M 

s 

•  The  ratio  of  maximum  P  to  P^  amplitudes  is 
a  function  of  distance,  with  maxima  generally 
falling  in  the  range  of  300  to  750  km,  although 
considerable  variation  is  evident  from  station  to 
station.  These  observed  differences  lead  to  the 
conclusion  that  m^  is  currently  unreliable  when 
based  on  data  recorded  less  than  1000  km  from 
the  source 
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•  Nonlinear  or  "exact"  confidence  regions  may  be 
computed  from  use  of  the  last  equations  in  Section  II, 
subsection  B 

•  Chengmai,  Thailand,  is  the  most  capable  single-element 

station  studied,  having  a  uieoretical  capability  to  record 

all  events  of  m,  >  4.  5  to  more  than  100°.  Another 
b  — 

1 4  stations  indicate  capabilities  to  record  events 

from  m.  =  4.  5  to  5.  0  at  similar  distances:  five  of 
b 

these  are  Arctic  or  Antarctic  stations;  three  are 
located  in  the  Western  United  States;  three  are  on 
the  Indian  subcontinent;  two  are  in  Africa;  and  the 
ether  is  in  Shiraz,  Iran 

•  /Analysis  of  magnitude  residuals  indicates  that 
patterns  of  residuals  do  exist;  these  patterns  might 
be  used  to  infer  source  mechanisms  and,  possibly, 
as  a  further  criterion  for  discriminating  between 
earthquakes  and  explosions 

Results  included  in  the  appendixes  can  be  summarized 

as  follows: 

•  Approximately  half  of  the  344  January  1964  revised 
hypocenters  are  changed  appreciably  from  the 
USC&GS  preliminary  locations 

•  Depth  determination  remains  a  considerable  prob¬ 
lem  and  requires  continued  investigation 
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•  The  number  of  earthquakes  occurring  in  the 
vicinity  of  Istanbul,  Turkey,  increases  much 
more  .lowly  with  decreasing  magnitude  than 
the  number  near  Shiiaz,  Iran;  Quetta,  Pakistan; 
or  Shillong,  India.  However,  more  larger 
earthquakes  occurred  near  Istanbul  than  near 
the  other  stations  in  the  time  period  studied 

•  Nearly  3-1/2  times  as  many  pP  phases  are 
identified  after  application  of  array  processing 
than  can  be  identified  from  the  single- channel 
presentation,  based  on  the  sample  investigated 

•  Application  of  digital  filters  provides  better 
signal  enhancement  than  does  analog  filtering 

•  Studies  of  energy  attenuation  based  on  automated 
calculation  of  energy  from  digital  data  indicate 
considerable  error  in  the  Q  curves  used  in  mag¬ 
nitude  calculation 

•  Empirical  fits  to  observed  energy  attenuation  as 
a  function  of  distance  lead  to  a  relationship  of  the 
form 

rp  v  -2  .  -2K  r 

E  =  K  r  cos  (e  )  e  2 
1  o' 

where  and  are  constants,  r  is  epicentrai 
distance,  and  e^  is  angle  of  incidence 

•  Assuming  that  the  shape  of  power  spectra  is 
independent  of  distance,  as  shown  by  LeBlanc 
and  Howell  (1965),  a  set  of  crustal  filters  is  de¬ 
veloped  to  describe  in  the  frequency  domain  the 
effect  of  crustal  structure  beneath  several  stations 
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These  results,  cor  ibined  with  those  previously  published, 
lead  to  the  following  conclusions. 

•  A  hypocentei  program  capable  of  high  accuracy 
should  include  the  means  of  correcting  travel-time 
tables  on  the  basis  of  both  station  and  source  legion 

•  High-quality  well-distributed  stations  yield  better 
results  than  do  large  numbers  of  average  stations 
concentrated  in  one  or  two  quadrants 

•  Generally,  magnitudes  computed  at  distances  of 
less  than  1000  km  are  inaccurate  and  should  not 
be  used  as  estimates  of  event  size  or  be  combined 
with  teleseismic  data  to  obtain  average  values 

•  Short-period  P-wave  magnitudes  now  in  use  can 
be  considered  accurate  to  only  bout  ±  0.  5  units; 
attempts  to  compare  m^  with  other  magnitudes 
must  consider  this  scatter 

•  Variations  in  seismic  activity  represent  more 
than  simple  differences  in  level;  the  ratio  of  the 
numbers  of  large  and  small  earthquakes  also 
varies  considerably  from  place  to  place 

Considering  the  results  and  conclusions,  the  following  recom¬ 
mendations  are  made: 

•  More  studies  of  crustal  structure  are  needed  in 
areas  of  prevalent  seismic  activity  as  well  as 
near  seismograph  stations;  such  studies  will 
benefit  both  hypocenter  and  magnitude  determina¬ 
tions 


1-4 


•done*  •orvioos  division 


•  There  should  be  more  use  of  tape  recording  — 
preferably  direct  digital  recording 

•  The  "world-array"  processing  concept  should 
be  applied  to  more  data  and  to  an  even  broader 
range  of  problems  than  those  discussed  in  this 
report 

•  Automated  techniques  for  "energy"  or  magnitude 
determination  should  be  considered  as  alternates 
to  current  methods 

j  The  method  for  obtaining  first  approximations  of 
hypocenter  coordinates  developed  under  th>  ’  con¬ 
tract  should  be  investigated  more  fully 

•  The  method  f^r  assessing  epicenter  accuracy 
(described  in  Special  Report  No.  5)  should  be 
developed  and  applied  more  fully 

•  The  USC&GS  should  provide  continuing  year-to- 
year  assessments  of  seismic  activity  on  regional 
and  worldwide  bases 

•  The  USC&GS  should  develop  as  complete  a  file 
of  seismic  data  as  possible  in  a  format  adopted 

as  "standard"  on  at  least  a  national,  and  preferably 
on  an  international,  basis 

These  recommendations,  if  followed,  would  provide  consid¬ 
erable  steps  toward  solving  several  problems  in  earthquake  seismology. 
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SECTION  II 

PRESENTATION  OF  RESULTS 


A.  MAGNITUDE  STUDIES 
1.  Magnitude  Calculation 
a.  Definitions 

Various  methods  for  computing  seismic  event  magnitudes  have 
been  developed  over  the  years,  beginning  with  Richter's  (1935)  local  magni¬ 
tude  scale  (M^).  Gutenberg  M5,  a,  b,  c)  extended  magnitude  determination 
to  20-sec-period  surface  waves  and  the  bod'  phases  P,  PP,  and  S.  Bath 
(1952)  extended  magnitude  computation  to  vertical-component  surface  waves 
and  broadened  the  period  range  over  which  surface-wave  calculations  could 
be  made.  Gutenberg  and  Richter  later  (1956)  revised  bodywave  magnitude 
calculation  and  the  relationship  between  surface-wave  magnitudes  (M  )  and 
what  Gutenberg  termed  "unified  magnitude"  (m). 

The  USC&GS  (1963)  adopted  a  form  of  Gutenberg's  unified 
magnitude  (m^)  and  began  reporting  values  of  this  magnitude  in  its  Prelimi¬ 
nary  Determination  of  Epicenter  cards  in  April  1963.  This  latter  magnitude 
(m^)  is  calculated  by 

mb  =  log  Y  +  °  (1) 


where 


A  =  maximum  amplitude  of  Pn  or  P  in  the  first  few 
cycles,  measured  in  microns,  as  recorded  by 
short-period  vertical- component  seismographs 

T  =  dominant  period  of  the  measured  P  wave 

Q  =  depth-distance  factor  as  given  by  Gutenberg  and 
Richter  { 1956)  for  the  distance  range  5°  to  110° 
and  extrapolated  to  2°  for  surface-focus  events, 
assuming  an  inverse  cube  attenuation  of  signal 
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Texas  Instruments  Incorporated  has  used  all  of  the  preceding 
methods  of  magnitude  calculation.  Taking  into  account  previously  observed 
variations  in  similar  magnitudes  computed  from  data  recorded  by  long-period 
and  short-period  instruments,  the  following  10  magnitudes  were  computed, 
here  possible,  from  1964  data: 


M 


1 


M 

8 

M9 


M 


10 


short-period  vertical-component  P-wave 
magnitude  differing  from  only  at  short 
distances  (A^  10°)  where  the  amplitude 
used  is  the  maximum  P  amplitude  (not 
necessarily  P  ) 

long-period  vertical-component  P-phase 
magnitude 

short-period  vertical-component  PP-phase 
magnitude 

long-period  vertical- component  PP- phase 
magnitude 

short-period  horizontal- component  S-phase 
magnitude 

long-period  horizontal- component  S-phase 
magnitude 


vertical- component  surface-wave  magnitude 
M 

s 


b.  Method 


The  first  seven  magnitudes  are  computed  from  bodywave  data 
computed  from  equations  similar  to  (1): 


M. 


A. 

log  y~  + 


i  =  1-7  (2) 


n-2 
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where 


vertical- component  amplitude  (peak-trough/2) 
of  P  or  PP  phase,  corrected  for  instrument 
response,  divided  by  the  dominant  period  (|i/sec) 
for  i  =  1  to  5;  for  i  =  6  and  7,  the  vector  ampli¬ 
tude  of  SH  is  used 

depth-distance  factor;  for  i  =  1  to  3,  PZ  values 
are  used;  for  i  =  4  and  5,  PPZ  values  are  used; 
for  i  =  6  and  Z,  SH  values  are  used;  all  values 
are  as  given  by  Gutenberg  and  Richter  (1956) 


The  eighth  magnitude  is  an  adaptation  of  Richter's  (1935)  local 

magnitude  scale  M  .  It  is  computed  by 
L 


M8  =  iog  \A 


-  log  A 


where 


maximum  short-period  horizontal  trace  motion 


■  log  = 


response  of  standard  Wood-Anderson  seismograph 
at  period  T 

response  of  short-period  seismograph  recording 
the  data  at  period  T 

Richter's  (1935)  distance  factor  extrapolated 
from  600  km  to  10° 


Separate  magnitudes  are  computed  for  each  horizontal  com¬ 


ponent  and  the  mean  taken. 


The  ninth  and  tenth  magnitudes  are  computed  from  surface- 
wave  data  in  the  period  range  of  17  sec  <  T  <  23  sec,  using  the  following 
equations,  adapted  from  B£th(1952)  and  Gutenberg  (1945a): 

M9  =  a  [log  A  +  -|-  log  +  24.  13  A  (KT~  K^) 

+  b  -  log  bJ  +  C 

and  (4) 

1  20 

M1q  =  log  A  +  j  log  —  +  24.  13  A°  (KT-  K^)  -  log  B 


where 

A  =  maximum  amplitude  of  surface  waves 

(vector  sum  of  both  horizontal-component 
amplitudes  for  MjQ)  in  n 

T  =  period  of  surface  waves  measured 

(17  sec  S  T  s  23  sec)  -  log  B  =  distance 
factor 

A  =  epicenter-to-station  distance  in  degrees 

(K^,  -  K^q)  =  difference  in  extinction  coefficients  of 
surface  waves  of  T-sec  and  20-sec 
periods 

h  =  depth  of  focus  in  km 

a  =  0.8  for  h  5  40  km,  0.  5  for  40  km 

<  h  s  300  km 

b  =  0.  0082  h  for  h  £  100  km,  0.85  for 

100  km  <  h  5  300  km 

c  =  1.42  for  h  s  40  km,  3.55  for  40  km 

<  h  5  300  km 
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Values  of  are  calculated  only  for  events  occurring  at 
depths  £  35  km  and  for  events  occurring  at  depths  ^  300  km. 

A  computer  program  described  m  Semiannual  Technical 

Report  No.  Ill  (Texas  Instruments,  1966)  was  written  to  compute  earthquake 

magnitudes.  This  program  computes  each  of  the  10  magnitudes  for  each 

station  recording  a  particular  event  (where  data  permit).  The  mean  and 

standard  deviation  of  each  magnitude  type  for  each  event  also  is  computed 

where  data  permit.  Values  of  event  mean  magnitudes  (m,  =  M,  and  M  = 

bis 

M^)  are  given  in  Appendix  A  with  results  of  hypocenter  determinations. 

2.  Comparison  of  Magnitude  Scales 
a.  General 

In  the  previously  cited  technical  report  by  Texas  Instruments, 
various  magnitudes  were  compared.  From  the  standpoint  of  frequency  of 
calculation,  it  was  determined  that  short-period  P-wave  magnitudes  were 
most  valuable  since,  in  the  study  of  1963  seismicity  (Texas  Instruments, 
1964b),  more  than  2-1/2  times  as  many  short-period  P-wave  magnitudes 
were  computed  than  all  other  types  combined.  The  internal  variation  was 
greatest,  however,  for  the  short-period  P-phase  magnitudes  than  for  all 
others  except  short-period  S-phase  magnitudes. 


Variations  in  the  short-period  P-phase  magnitude  could  be 
reduced  perhaps  by  revising  the  depth-distance  factors  (Q  tables).  However, 
it  has  been  observed  that  the  m^  magnitudes  do  not  appear  to  behave  properly 
for  large  events;  furthermore,  data  for  computing  m^  for  rather  large 
earthquakes  are  difficult  to  obtain  since  the  dynamic  range  of  the  seismo¬ 
graphs  is  often  exceeded,  precluding  measurement  of  period  and  amplitude. 
To  alleviate  this  situation  and  to  provide  consistent  magnitude  data  for  both 
the  smallest  and  largest  events  recorded,  it  is  desirable  to  use  an  alternate 

magnitude  scale  which  may  be  converted  to  m  values. 

b 
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b.  m  and  M  (M  and  M.  .) 
b  si  1 0 

In  Technical  Report  No.  Ill  (1966),  the  better  1963  magnitude 

data  were  investigated  to  obtain  a  relationship  between  m  and  M  .  Actually, 

b  s 

two  relationships  were  found.  When  m  was  considered  independent,  the 

b 

relationship  was  determined  to  be 


M  =  0. 93  rrv  +  0.565 
s  b 


When  M  was  considered  independent,  the  derived  relationship 


rm  =  0.439  M  +  2.825 
b  s 


Neither  relationship  was  in  good  agreement  with  Richter’s 
(1958)  relationship: 

M  =  1.59  m  -  3. 96 
s  b 


m,  =  0. 63  M  +2.5 
b  s 


Visual  inspection  of  the  data  used  to  determine  the  m,  -  M 

b  s 

relationship  indicated  that  a  line  with  the  same  slope  as  the  Richter  relation¬ 
ship  appeared  to  fit  the  data  as  well  as  the  relationships  given  by  Equations 
(5)  and  (6).  Consequently,  a  new  relationship  was  derived,  assuming  a 
slope  of  1.59.  Resulting  from  this  was  the  relationship 

M  s  1.59  m  -  2.82 
s  b 


m,  =  0.63  M  +  1.77 
b  s 
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Figure  II-  I  shows  the  three  derived  relationships  and  the  data 
used  in  their  derivation. 

c.  m^  and  Long-Period  P- Phase  Magnitudes  (M^  and  '"  ) 

Linear  reg  ession  of  (long-period  P-wave  magnitude)  on 

M  (m  ),  using  1963  data,  obtained  the  relationship 
1  b 

M3  =  0.83  +  1.67  (9) 

Since  the  slope  was  nearly  unity,  the  differences  between  M„ 
and  ir»k  were  investigated  to  determine  average  differences  for  each  seismo¬ 
graph  station  and  the  mean  difference  for  the  total  data.  Table  II- 1  lists  the 
corrections  to  be  applied  to  values  for  each  station  to  obtain  an  estimated 
m^  and  the  mean  correction  for  all  stations.  Using  this  mean  value,  may 

be  related  to  m_  by 
b 

M  =  m  +  C.  72  (10) 

3  b 

Figure  II- 2  shows  the  relationships  given  by  Equations  (9)  and 
(10)  and  the  data  used  in  their  derivation. 

It  is  doubtful  that  significant  differences  in  internal  consistency 
would  be  noted  in  converting  values  to  the  m>  scale  by  any  of  the  three 
following  methods: 

•  To  average  values  computed  at  all  stations 

and  to  convert  to  m  using  Equation  (8) 

b 

•  To  average  all  values  for  the  event  and  to 

convert  to  m  using  Equation  (9) 
b 

•  To  apply  the  corrections  listed  in  Table  II- 1  to  the 

values  computed  from  each  station's  data  and 
then  to  average 
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. . 


As  shown  in  Table  II-2,  however,  there  appears  to  be  some 
correlation  between  geographical  location  and  the  difference  between  and 
m^.  Thus,  the  third  approach  may  be  the  best  method  for  approximating  m^ 
from  long-period  P-wave  magnitudes. 

3.  Horizontal-  and  Vertical-Component  Surface-Wave  Magnitudes 
(M1q  and  Mg) 

Vertical -component  surface-wave  magnitudes  generally  are 
higher  than  horizontal- component  surface- wave  magnitudes.  This  may  result 
from  deletion  of  Bath’s  (1952)  path  correction  in  Equation  (4).  The  correction 
was  assumed  to  be  0,  since  no  prior  information  was  available  for  the  stations 
used. 


In  Texas  Instruments  Technical  Report  No.  Ill,  the  results  of 
a  linear  regression  of  Mg  on  M^  were  reported  as 


Mg  =  0.764  M1q  +  1.728 

However,  the  Mg  scale  was  originally  defined  in  such  a  way 

to  make  it  equal  to  the  M  scale.  Therefore,  the  differences  between  M 

s  9 

and  M  calculations  have  been  tabulated,  and  mean  corrections  of  M  to 
s  9 

obtain  M  values  have  been  determined.  The  results  are  shown  in  Table  II- 1. 
s 

The  mean  difference  of  all  and  M  values  also  has  been 

9  s 

determined.  However,  a  comparison  a  the  standard  deviations  for  the  data 
of  individual  stations  and  for  the  total  mean  correction  indicates  that  the 
individual  corrections  are  probably  preferable.  Also,  Figure  II- 3  indicates 
considerable  geographical  dependence  for  the  station  corrections. 


j 


I 
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Magnitude  Deviation  from 


4.  Comparison  of  P  and  Maximum  P  Amplitudes,  A  <  9° 

n 

As  discussed  earlier,  the  scale  is  a  somewhat  specialized 
form  of  Gutenberg  and  Richter*s  {1956)  unified  magnitude  scale.  However, 
the  definition  of  the  USC&tGS  differs  in  one  respect  from  the  Gutenberg  and 
Richter  P-phase  magnitude:  P^  amplitudes  rather  than  maximum  P  amplitudes 
are  used. 

At  distances  of  less  than  about  1000  km  or  9°,  considerable 
variation  in  magnitudes  computed  by  the  two  methods  has  been  observed. 

Since  the  two  magnitudes  differ  only  in  the  amplitude  used,  comparison  results 
in  an  analysis  of  P^  and  maximum  P-amplitude  variations. 

Using  data  accumulated  in  the  1963  study  (Texas  Instruments, 
1964b),  variations  between  M^  and  M^  values  for  the  following  stations  have 
been  analyzed: 

•  Istanbul,  Turkey  (1ST) 

•  Shiraz,  Iran  (SHI) 

•  Quetta,  West  Pakistan  (QUE) 

•  Shillong,  India  (SHL) 

•  Chengmai,  Thailand  (CHG) 

•  Port  Moresby,  New  Guinea  (PMG) 

•  Nana,  Peru  (NNA) 

•  Arequipa,  Peru  (ARE) 

•  Antofagasta,  Chile  (ANT) 

Both  M  and  M  magnitudes  were  computed  a  total  of  1566 

L  w 

times  at  all  stations.  About  50  percent  of  the  time,  the  analysts  measured 
the  same  amplitude  and  period  for  both  magnitude  calculations.  At  least 
two  interpretations  of  this  fact  may  be  made: 

•  About  half  the  time,  Pn  is  the  maximum 
P  phase  at  distances  ^  1000  km 

•  Often,  the  Pn  phase  is  not  measurable 


11-14 


selsnes  ssrvless  division 


Th  .  experience  of  Texas  Instruments  in  analyzing  earthquake 

and  explosion  data  over  the  past  6  years  indicates  that,  at  most  distances  for 

most  stations,  P  is  not  the  maximum  P  amplitude.  Thus,  the  second  inter- 
n 

pretation  is  preferred. 

Of  the  62  events  within  luCO  km  of  Istanbul  which  were  analyzed, 
three  had  epicenters  located  by  the  USC&GS;  for  the  remainder,  only  the  epi- 
central  distance  determined  from  the  S-P  interval  was  known.  and 

magnitudes  were  equal  for  30  of  the  62  events  (48  percent).  Figure  II-4  shows 
the  differences  in  magnitudes  for  the  other  32  events.  At  distances  of  approxi¬ 
mately  1°  to  2°,  maximum  P  amplitudes  averaged  about  three  times  the  P 

n 

amplitudes.  From  approximately  2°  to  3°,  maximum  P  amplitudes  were 
about  twice  as  large  as  P^  (on  the  average).  Little  data  are  available  beyond 
3°  distance,  but  there  are  indications  that  the  ratio  of  maximum  P-to-P 

n 

amplitudes  increased  to  approximately  5  around  a  distance  of  5°.  At  greater 
distances,  the  ratio  approaches  unity. 

At  Shiraz,  Iran,  only  25  of  the  99  analyzed  events  had  = 

The  majority  of  the  other  74  events  were  at  1  °  to  2°  from  the  station,  and 
considerable  scatter  was  noted  from  these  data.  However,  an  average  ratio 
of  approximately  3  to  4  is  indicated  for  the  distance  range  of  1  0  to  2°.  The 
ratio  appears  to  increase  to  approximately  10  near  4°  and  slowly  decrease 
to  near  0  at  9°  (1000  km),  as  shown  in  Figure  II- 5. 

The  two  magnitudes  were  equal  for  approximately  one-third 
of  the  analyzed  events  near  Quetta.  The  computed  ratios  for  the  remainder 
of  these  events  exhibited  a  distinctly  different  pattern  from  the  previous  two, 
as  shown  in  Figure  II -6 .  In  this  case,  maximum  P  amplitudes  average  about 
three  times  the  P^  amplitudes  from  approximately  1°  to  3°.  The  ratio  de¬ 
creases  to  about  2  at  a  distance  of  4°  and  approaches  0  at  approximately  5°. 

At  about  6.5°,  however,  the  average  ratio  sharply  jumps  to  a  value  of 
approximately  9  and  decreases  to  about  2  again  near  7.5°.  Little  data  are 
available  at  greater  distances,  but  the  ratio  does  not  seem  to  approach  1.0  at 
9°  as  was  exper.ed. 
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gure  II- 5.  Comparison  of  Maximum  P  and  Pn  Amplitudes 
in  Vicinity  of  Shiraz,  Iran 
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Insufficient  data  were  obtained  for  valid  analyses  of  P  and 

n 

maximum  P  amplitudes  at  Shillong,  India,  and  at  Chengmai,  Thailand,  as 
two-thirds  of  the  events  at  Shillong  and  half  of  the  events  at  Chengmai  had 
equal  magnitudes  on  each  scale.  From  the  little  data  available,  however, 
the  average  amplitude  ratio  for  both  stations  can  be  estimated  to  be  approxi¬ 
mately  2  for  most  of  the  distance  range. 

CVer  40C  events  at  Port  Moresby,  New  Guinea,  were  analyzed 

in  this  study.  Since  the  noise  level  at  this  station  is  high,  it  is  not  surprising 

to  find  that  nearly  half  of  the  events  had  equal  magnitudes  on  both  scales. 

Considerable  scatter  was  noted  over  the  entire  distance  range,  even  for 

events  with  known  epicenters.  This  scatter  complicated  the  determination 

of  any  clear  pattern  of  amplitude  ratios.  High  ratios  appeared  prevalent  from 

approximately  3°  to  5°  (maximum  P  up  to  12  times  P  amplitudes).  From 

n 

about  6°  to  8°,  most  events  were  located  by  the  USC&GS;  the  sverage  ratio 
over  this  range  was  approximately  4.  Relatively  high  amplitude  ratios  were 
observed  to  9°.  Figure  II- 7  shows  the  distribution  of  amplitude  ratios  with 
distance. 

Slightly  over  40  percent  of  the  analyzed  events  near  Nana, 

Peru,  had  equal  magnitudes  on  both  scales.  The  pattern  here  seems  to 
have  been  increasing  ratios  to  an  average  value  of  about  4  at  a  distance  of  3° 
and  near  this  level  to  approximately  7°,  as  shown  in  Figure  II-8.  Data 
were  sparse  at  the  greater  distances. 

More  than  half  of  the  500  analyzed  events  near  Arequipa, 

Peru,  had  equal  magnitudes  on  both  scales.  The  average  amplitude  ratio 
was  2  to  3  at  1°  to  7°,  with  an  indication  of  an  approach  to  unity  near  9°* 
as  indicated  in  Figure  II- 9 . 
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Figure  II- 7.  Comparison  of  Maximum  P  and  P  Amplitude 
Vicinity  of  New  Guinea 


Comparison  .1  Maximum  P  and  Pn  Amplitudes 
in  Vicinity  of  Nana,  Peru 
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in  Vicinity  of  Arequipa,  Peru 


Data  near  Antofagasta,  Chile,  were  from  events  in  a  narrow 

distance  range  of  approximately  1.5°  to  3"'  The  average  amplitude  ratio 

was  about  2.  5  to  3  over  this  distance  interval.  Of  the  events  analyzed  in  this 

area,  75  percent  had  equal  magnitudes  on  both  scales.  This  probably  resulted 

from  high  noise  levels  and  low  instrument  gain  which  'creases  the  probability 

of  P  detection, 
n 

In  summary,  maximum  P  amplitudes  at  most  stations  appear 
to  be  at  least  twice  as  large  as  P  amplitudes.  High  amplitude  ratios  often 
occur  from  3°  to  5C  distances.  Variations  in  patterns  rom  station  to  station 
indicate  that  pa  rameters  other  than  distance  affect  propagation  patterns  for 
various  branches  of  P.  Some  of  these  parameters  might  include  azimuth, 
source  region,  crustal  structure,  travel-time  curves,  and  source  mechanism. 
Investigation  of  the  effects  of  these  parameters  is  needed  before  reliable 
magnitudes  can  be  computed  at  sh.Trt  distances. 

B.  EXACT  CONFIDENCE  REGIONS 

1.  Definitions  ?nc’  A  ssumptions 

a.  Propagation  of  Seismic  Waves 

For  che  purposes  of  mathematical  characterization,  the  physical, 
or  spatial,  earth  is  defined  as  follows: 

•  The  spatial  earth,  denoted  by  V,  is  a  compact, 
convex  continuum  with  the  usual  3 -dimensional 
metric  and  has  the  geometric  and  topological 
properties  of  the  closure  of  the  interior  of  an 
oblate  spheroid  except  that  there  is  a  positive 
real  number  t  such  that,  for  any  boundary  point 
p,  the  distance  from  p  to  the  center  of  W  may 
differ  from  the  analytic  value  by  an  amount,  the 
absolute  value  of  which  is  less  than  or  equal  to  t 

•  The  surface  of  the  spatial  earth,  denoted  by  B(W), 
is  the  boundary  set  of  W 
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Regarding  the  temporal  continuum  as  a  copy  of  the  real  line 
denoted  by  T,  the  spatio-temporal  earth  is  defined  as  follows: 

e  The  spatio-temporal  earth  T  x  W  is  the 
4-dimensional  space  generated  by  T  and 
directed  by  W 

Under  these  conditions,  a  seismic  event  occurring  at  time  H 

o 

and  having  the  spatial  (curvilinear)  coordinates  X  ,  $  ,  and  R  can  be  treated 

o  o  o 

just  as  a  vector  (H,X,$,R)eTxW.  Now,  let  (u>  ,  ....  u>  )  be  a  finite 
o  o  o  o  1  y 

set  of  wave  phases  associated  with  the  seismic  event  (H  ,  X  ,  $  ,  R  )  and 

[o  o  o  o 

H  ,  H  '  ,  denoted  by  T  ,  be  a  closed  interval  on  T.  Then,  the  basic 
.  o  o  o 

assumption  concerning  the  propagation  of  seismic  waves  can  be  stated  at 
follows: 

(1)  For  each  phase  uu^,  there  is  an  associated  continuous 

monotone  transformation  0. :  W  xl(H  ,  ,  ,  R  )f  -  T 

1  '  o  o  o  o'  o 

so  that  a  point  fx...  ♦  ,  R,  (H  ,  X  ,  $  ,R  j|  e  Cl. - 1  (H)  where 
L  o  o  o  oj  i 

H  S  Hs  H  if,  and  only  if,  H  is  the  time  at  which  the 
o  o 

<u.  phase  arrives  at  (X,  I,  R) 

In  accordance  with  the  characterization  given  in  assumption  (1), 

the  point  (H  ,  X  ,  ?  ,  R  )  which  occurs  in  the  definition  of  the  domain  of  Q 
o  o  o  o  i 

is  referred  to  as  the  parameter  of  the  transformations.  Alternately,  equivalent 

transformations  could  have  been  developed  as  temporal  metric  functions  on  the 

space  T  x  W;  under  such  functions,  the  distance  between  two  points  (H  ,  X  , 

o  o 

$  ,  Rq)  and  (H,  X,  $,R)  would  simply  be  the  absolute  difference  between  origin 
time  Hq  and  the  U).  phase  arrival  time  H  at  (X,  i,  R).  However,  the  approach 
taken  in  the  statement  of  the  assumption,  with  the  consequent  representation 
of  T  x  W  as  a  parameter  space,  is  better  suited  to  the  statistical  development 
which  is  to  follow. 
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It  should  be  noted  that  assumption  (I)  does  not  strictly  agree 
with  empirical  evidence,  because  travel-time  curves  for  certain  wave  phaseo 
show  some  signs  of  looping.  Hence,  for  regii  na  in  the  domain  corresponding 
to  a  loop  in  the  graph  of  fX,  the  assumptions  that  Q.  is  single- valued  and 
monotone  are  unwarranted.  In  practice,  this  difficulty  can  probably  be  over¬ 
come  in  the  case  of  large  samples  by  existing  routines  for  elimination  of  data 
differing  too  much  from  the  sample  mean;  in  the  case  of  small  or  minimal 
samples,  it  may  be  necessary  to  include  a  procedure  for  eliminating  obser¬ 
vation  points  lying  in  the  region  of  a  loop  as  part  of  the  process  of  refinement. 

b.  Analytic  Description 


The  combination  of  tables  giving  wave-phase  travel  times, 

depth  corrections,  station  corrections,  and  ellipticity  corrections,  together 

with  rules  for  using  these  tables,  can  be  thought  of  as  defining,  for  each  phase 

;ju.,  a  function  f. :  W  x  (T  x  W)  -  T  so  that  f.  [>,,  J,  R,  (H  ,  X  ,  f  ,  R  )1  is  the  com- 
i  1  iL  ooooJ 

puted  time  at  which  phase  uu  should  arrive  at  (X,  i,  R).  Taking  this  view,  the 
basic  assumption  regarding  the  analytic  description  of  the  propagation  of 
seismic  waves  can  be  stated  as  follows: 


(2) 


<(X.,f.,Rj):U  j 
points  (stations); 


£n|  is  an  arbitrary  set  of  observation 
therefore,  .e.,  the  error  residual  at 


the  j  point  with  respect  to  wave  phase  u'.,  is  given  by 

the  statement  .  e.  =  Q.[\.,$.,R.,  (H  ,  X  ,  £  ,R  )"| - 
ij  i  jj  j  ooo  o  J 

f.  [X R.,  (H  ,X  ,  *  ,R  )"];  the  error  residuals  .  e. 
i  jjjoooo-1  1J 

are  independent  normal  variates  with  mean  zero 

2 

and  variance  a 


The  importance  of  assumption  (2)  can  be  discussed  more  mean¬ 
ingfully  in  the  light  of  the  statistical  theory  which  is  to  be  developed;  accordingly, 
the  discussion  is  postponed  to  the  conclusion  of  that  development.  It  should  be 
understood  from  the  outset,  however,  that  this  is  the  central  assumption.  The 
justificatior  of  the  present  approach  turns  on  the  validity  of  this  assumption. 
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c.  Hartley's  Modified  Gauss -Newton  Method 

The  procedure  for  nonlinear  hypocenter  determination  to  be 
described  in  the  sequel  employs  Hartley's  modification  of  the  Gauss-Newton 
nonlinear  estimation  method.  In  addition  to  assumption  (2),  Hartley's 
method  proceeds  from  the  following  definitions  and  assumptions.  (For 
clarity,  subscripts  which  serve  only  to  associate  an  expression  with  a  partic¬ 
ular  wave  phase  will  be  used  from  this,  point  forward  only  when  the  distinction 
is  critical. )  Let 

C  5  (C1’C2’C3-C4) 

where 


let 


where 


For  a  given  set  j(X.^,  RJ  :  1  <  i  <  nj  of  observation  points, 


X.  =  R. 
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(3)  For  each  i,  jeF,  £.  (X,  C)  and  f. .  (X,  C)  are  continuous 

1  ij 

functions  of  C.  and  C  for  all  X,  where  keN 
i  J  k 

(4)  For  any  rontrivial  set  of  real  numbers  ju.  lieFlwith 

2  1 
Z.  „(u.)  >  0,  it  is  the  case  that 

ie  F  i 

Z,  KT  (Z.  _  [u.f.  (X,  ,  C)1  j  >  0 
keN  \  ieF  L  i  i  k  J  j 

for  all  C  lying  in  a  bounded  convex  set  S  of  the  param¬ 
eter  space  T  x  W 
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(5)  Let~S  be  the  complement  of  the  set  S  described 
in  assumption  (4)  and  let 


Q  =  lim  inf  Q  (X,  C) 

~  S 

Then,  there  exists  a  point  qC  lying  in  the  interior 
of  S  so  that 

Q  (X,  C)  <  Q 
o 


As  with  assumption  (2),  the  discussion  of  the  justification  for 
assumptions  (3),  (4),  and  (5)  is  best  deferred  to  the  conclusion  of  the  descrip¬ 
tion  of  the  estimation  technique.  The  purposes  served  by  these  latter  assump¬ 
tions  are  fairly  straightforward.  Assumption  (3)  makes  possible  the  def¬ 
inition  of  the  first  and  second  partial  derivatives  of  function  Q  with  respect  to 
Q  and  C.  where  i,  jeF;  these  definitions  will  be  stated  explicitly  in  the  sequel. 
According  to  Hartley,  assumption  (4)  is  "equivalent  to  the  well-known  assump¬ 
tion  of  nondegeneracy  of  rank  in  a  linear  least-square  problem.  .  .  "  and  is 
"usually...  satisfied  in  practical  situations."  The  modified  Gauss-Newton 
method  is  an  iterative  process,  and  the  role  of  assumption  (5)  is  to  guarantee 

convergence  of  the  iteration,  given  a  vector  C  of  the  specified  sort  as  a 

o 

starting  point. 

As  will  be  seen  later,  assumption  (4)  must  be  strengthened 
to  insure  that  the  solution  to  which  the  process  converges  (say,  C#)  gives  a 
unique,  absolute  minimum  for  Q(X,C). 

2.  Procedure 

The  following  conventions  are  introduced  so  that  the  development 
might  be  expressed  in  terms  of  matrix  algebra. 
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(1)  Let  i,  jeF,  not  necessarily  distinct,  and  let  kcN. 
Then, 

•  The  symbol  fl(X,C)  denotes  the 
N  x  1  column  vector  having 

f.  (X  ,  C)  as  its  It**1  component 

i  K 

•  The  symbol  f  (X,  C)  denotes  the 

N  x  4  column  vector  having  the 

th 

element  f.  (X,  ,  C)  in  the  k  row 

i  k 

and  the  j*"  column 

(2)  The  symbol  Y  denotes  the  N  x  1  column  vector  with 

th 

the  k  component  Y  ~  Q(X  ,  C) 

K  K 


Under  definitions  (1)  and  (2), 

Q  (X,  C)  =  [Y  -  f  (X,  C)]1  [Y  -  f  (X,  C)] 

=  Y*  Y  -  2Yt  f  (X,  C)  +  f  (X,  C)1  f  (X,  C) 

Hence,  for  i€ F, 


Q'tX.C)  sjc-  Q  (X,  C)  = 
i 


-2YV(X,C)  +  f(X,C)V  (X,  C)  4  fl  {X,C)t{(X,C) 
-2yV(X,C)  +  2f(X,C)V(X,C) 


=  - 2[y 1  -  f(x,  c) t]  ^(x,  c; 


=  - 2  [Y  -  f(X,C)]tf1(X,  C) 


Now,  let  the  symbol  denote  the  least-squares  estimate  of 
the  hypocenter  location  determined  by  the  existing  program  under  the  assump¬ 
tion  of  linearity. 
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By  Taylor's  theorem,  f(X,  C  +  D  )  can  be  approximated 

o  o 


f(X,  C)  +  f.  (X,  C)  D 
o  loo 


“.here  is  the  1x4  vector  having  the  i  component 


6  C. 

D°  =  ~7“^ 
1  6  v 


fo  0  £  v  £  1.  The  error  of  the  approximation  is  given  by 


4  ff.  (X,  C  +  v'D  )D  f 
2  L  1  o  o  oj 


for  some  v'  so  that  0  i  v'  s  1  where  the  exponent  indicates  the  usual  formal 
expansion  of  the  operator. 

Substitution  of  the  approximation  for  f  (X,  C)  yields 


Q1  (X,  C)  =  -  2|y  -  f(X,  C)-f.(X,  OD1]  f'tX,  C)  (13) 

o  L  o  1  o  oj  o 

=  -  z[y  -  f(X,  C)]1  fl(X,  C)  +  2  D  f  (X,  CjVfX,  C) 
L  oj  o  oio  o 

=  Q1  (X,  C)  +  2fl  (X ,  C)*  f.  (X,  C)  D  1 
o  o  1  o  o 

The  formal  least- square s  equations  which  minimize  Q(X,n) 
are  given  by  the  system  Ql  (X,  C)  =  0  {i  =  1,  .  .  . ,  4);  hence,  for  i  =  1 ,  .  . .  ,  4, 


2f1(X,  C )tf.(X,  C)  D  1  =  -  Q1  (X,  C)  (14) 

o  1  o  o  o 
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By  assumption  (4),  the  determinant  of  the  latter  system  is 
of  rank  4;  therefore,  it  is  possible  to  solve  the  system  for  the  vector  D^. 

Now,  let 

Q  (v*)  =  Q  (X,  C  +  v'D  ) 

a  a  a 


for  0  <  v1  £  1  where  a  is  a  nonnegative  integer  and  let  v^  denote  the  value  of 
v1  for  which  (v')  is  a  minimum  on  the  interval  0  <  v1  ^  1.  Hartley  gives 

vi  =  T  +  T  [Qa(0>  •  Qa<“)]/  [Qa(U  -  2Qa(!)  +Qa(0>]  (15) 


The  estimate  giv  n  by  the  a  iteration  of  Hartley's  method 


is  defined  to  be  the  vector 


.  ,C  a  C  +  v'  D 

a  +  1  a  a  a 


If,  after  m  steps,  Hartley's  method  yields  an  estimate  C*  which 
satisfies  predetermined  convergence  criteria,  an  exact  lOOp  percent  con¬ 
fidence  region  for  Osps  1  can  be  constructed  by  Booker's  method.  As 
before,  the  linear  terms  of  a  Taylor  expansion  of  f(X,C)  about  C*  are 
given  by  f  (X,  C*)  +  f^  (X,  C*)  D*.  Let 


A  =  f  t  (X,  C*)  [fj  (X,  C*)  *  ^(X.C*)]’1  fL  (X.C*)1 


Reg(e)  =  [Y  -  f(X,C)]  A  [Y  -  f(X,C)] 


Res  (e)  =  [Y  -  f  (X,  C)]  (I  -  A)[Y  -  f(X,C)3 
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Then,  an  exact  iOOp  percent  confidence  region  is  given  by 

R  =  {CeT  x  W  :  Reg(e)/Res(e)  <  45(Pi4,N-4)/(N-4)f  (19) 

p 

where  Reg  (e)/Res  (e)  has  the  distribution  4$  (p;4,  N-4)/(N-4). 

C.  INSTRUMENTAL  PERCEPTIBILITY 

Instrumental  perceptibility,  a  measure  of  a  seismograph 
stations  capability  to  record  seismic  events,  is  based  on  noise  measurements 
and  theoretical  variations  in  P-phase  amplitudes  with  distance.  Before 
theoretical  perceptibility  Limits  can  be  calculated,  microseism  particle 
velocities  (amplitude /apparent  period)  for  each  station  must  be  estimated. 

This  is  accomplished  by  determining  the  distributions  of  seismic  noise 
(mM/sec)  at  each  station  and  using  as  the  noise  estimate  that  value  for  which 
the  probability  of  occurrence  is  0.  5.  Then,  the  observed  noise  at  a  particu¬ 
lar  station  should  be  less  than  or  equal  to  the  estimated  value  half  the  time. 

Solving  the  magnitude  equation  for  the  depth-distance  factor, 

one  obtains 

Q  =  mb  -  1o8(y)  =  f(A,h)  (20) 


Let  h  =  25  km,  substitute  estimated  values  of  A/T,  and  let 

=  4.  0,  4.  5  and  5.  0;  then,  the  distance  ranges  for  which  Q  £  the  value 

calculated  from  Equation  (20)  are  the  perceptibility  limits  for  m  =  4.  0, 

b 

4.5,  and  5.0  for  each  station.  The  results  of  these  calculations  for  44  of 


the  analyzed  stations  are  shown  in  Table  11-2.  These  results  indicate  that 
Chengmai,  Thailand,  has  the  highest  capability  for  recording  seismic  events 
but  is  followed  closely  by  two  Canadian  Arctic  stations,  Mould  Bay  (MBC) 
and  Resolute  Bay  (RES). 
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Table  II -2 

STATIONS'  CAPABILITIES  TO  RECORD  SEISMIC  EVENTS 


Station 

Y  |np/ sec) 

mb 

Instrumental  Perceptibinry  Limits  (°) 

4.  0 

0  -  2.7 

A  FI 

32 

4.  5 

0-4.0.  14.  3  -  10.  5 

5.  0 

0  -  6.0.  14.2  -  28.0.  39.8  -  12.5 

m 

0  -  2.  9 

AKU 

L.  1 

E9 

0  -  4.4,  15.2  -  20.0 

0  -  6.  3.  14.  1  -  29.6,  37.  2  -  43.9 

4.  0 

0  -  6.3,  14.0  -  30.  0,  36.8  -  44.  3 

ALE 

2.  5 

4.  5 

0-8.0,  12.  1-91.4 

5.  0 

0  -  104.2 

4.  0 

0  -  6.3,  14.0  -  30.0,  36.8  -  44.3 

ALQ 

2.  5 

4.  5 

0-8.0,  12.  1  -  91. 4 

5.  0 

0  -  104. 2 

0  -  4.  1,  15. 5  -  19. 7 

ARE 

9.5 

■9 

0  -  6.  1,  14.2  -  29.0,  39.0  -  43.3 

0-8.2,  12.  0  -  85.  0 

4.  0 

C  -  4. 1,  15. 3  -  19.  9 

ATL 

9 

4.  5 

0  -  6.2.  14.  1  -  29.4,  37.8  -  43.6 

5.  0 

0-8.5,  11.6  -  86.5 

4.  0 

0-2.7 

ATL 

35 

4.  5 

0-3.9,  15.8  -  19.  3 

5.  0 

0  -  5.8,  14.3  -  27.6 

4.0 

0-3.0 

BAG 

23 

4.  5 

0  -  4.5.  15.2  -  20.9 

5.  0 

0  -  6.5,  13.8  -  31.0,  36.4  -  45.0 

4.0 

0-4.5,  15.4  -  21.  1 

BOZ 

7 

4.  5 

0  -  6.2,  14.  1  -  29.  5,  37.  5  -  43.8 

5.0 

0  -  92.  0 

4.  0 

0  -  3.5 

BUL 

15 

4.  5 

0-5.1,  14. 6  -  23.  8 

5.  0 

0  -  7.  1,  13.  3  -  47.  8,  49.  5  -  56.  0,  71.2  -  81.8 

4.0 

0  -  7.  1,  13.3  -  47.8,  49.  5  -  56.0,  71.2  -  81.6 

CHG 

1.5 

4.5 

0  -  97. 0 

5.  0 

0  -  106. 0 
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Table  II- 2  (Contd) 
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Station 

Y  (tiiM / sec) 

mL 

Instrumental  Perceptibility  Limits  ( ° ) 

4.  0 

0  -  3.2 

KTG 

20 

4.  5 

0-4.7,  14.8  -  21.7 

5.  0 

0-6.7,  i  3.  7  -  32.  0,  35.  5  -  46.  0 

1 

0-4.7,  14.8  -  22.  0 

LON 

6 

H 

0  -  6.7,  13.6  -  33.  0,  35.  0  -  46.2 

0  -  93.  0 

1 

■ 

0-47,  14. 8  -  22. 0 

6 

0  -  '.7,  13.6  -  33.0,  35.0  -  46.2 

a 

0  -  93.  0 

0-2.2 

MAT 

66 

■9 

0-3.1 

‘■0 

0  -  4.6.  14.8  -  21.4 

4.  0 

0  -  6.7,  13.7  -  32.  0,  35.  5  -  46.  0 

MBC 

2 

4.  5 

0  -  92.  7 

5.  0 

0  105.2 

4.  0 

0-3.2 

MUN 

20 

4.  5 

0  -  4.7,  14  8  -  21.7 

5.  0 

0  -  6.7,  13.7  -  32.0,  35.5  -  46.0 

4.0 

0  -  5. 5,  14. 4  -  27. 0 

NAI 

4 

4.  5 

0  -  7.5.  13.  0  -  61. 0.  64.  0  -  83.  0 

5.0 

0  -  102.  0 

4.0 

0  -  5.C,  14.6  -  23.  5 

NDI 

5 

4.  5 

0  -  7.0.  13.  3  -  47.4,  50.  0  -  55.  5 

5.  0 

0  -  96.  5 

4.  0 

0  -  4.  5,  15.4-21.1 

NOR 

4.  5 

0  -  6.2.  14.  1  -  29.  5.  37.  5  -  43.8 

9.  0 

0  -  92.  0 

1.0 

0-4.0.  15.6  -  19.  5 

N'JR 

10 

4.  5 

0  -  6. 0,  14. 2  •  28. 0,  39. 5  -  43. 0 

5.  0 

0-8.0.  12.  1  -  84.  5 

4.  0 

0  -  4.  1,  15.  5  -  19.7 

PMG 

0.  T 

4.5 

0  -  6.  1,  14. 2  -  29.  0,  39.  0  -  43.  3 

PMG 


0  -  4.  1,  15.5-  19.7 

0-6.1,  14.4  •  £9.0,  19.0-43.3 

0-8.2,  12.  0  -  85.  0 


Table  II - Z  (Contd) 


Station 

Y  (mu/ see) 

mb 

Instrumental  Per.  bility  Limits  (  ) 

0-4.4,  15.0  -  20.6 

POO 

7.  5 

1 

0  -  6.4.  13.9  -  30.  5,  36.  6  -  44.  7 

1  1 

0-91.5 

0-5.6.  14.  3  -  27.  6 

QUE 

3.  5 

Bfl 

0-7.8,  12.  5  -  83.8 

5.  0 

0  -  102. 7 

mmm 

0  -  6.7,  13.7  -  32.  0.  35.  5  -  46.0 

RES 

2 

0-92.7 

r>i 

in 

o 

o 

0  -  2.9 

SCH 

25 

V i 

0-4.3.  15.0  -  20.3 

H 

0  -  6.  3,  14. 0  -  30.  0,  36.8  -  44.  3 

0-4.4.  15.0  -  20.6 

SCP 

7.5 

I  v-y- 

0  -  6.4.  13.9  -  30.  5,  36.6  -  44.7 

5.  o 

0-91.5 

4.  0 

0-5.8,  14.  3-27.6 

SHI 

3.  5 

4.  5 

0-7.8,  12.  5-83.8 

5.  0 

0  -  102.7 

4.  0 

0  -  6.  3,  14.  0  -  30.  0,  36.8  -  44.  3 

SHL 

2.  5 

4.  5 

0-8.0,  12.  1-91.4 

S.  0 

0  -  104.2 

I 

0-5.2.  14.  3  -  24.  1 

SPA 

4.  5 

IS  :  ^ 

0-7.2,  13.2-57.0,  71.0-  82.  1 

0  -  101.2 

D 

0-4.4,  15.0  -  20.6 

STU 

7.  3 

1 

0  -  6.4.  13.4  -  30.  5.  36.6  -  44.7 

1.  \ 

0  -  91.  5 

4.0 

0-3.5 

TRi 

in 

4.  5 

0  -  5.  1,  14.  6  -  23.  8 

5.0 

0  -  7.  1,  13.  3  -  47.8,  49.  5  -  56.  0,  71.2  -  81.8 

4.0 

0-5.0,  14.6-23.5 

WIN 

4.  5 

0  -  7.0,  13.  3  -  47.  4,  50  -  55.  5 

5.  0 

0  -  96.  5 
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Note  that  the  results  shown  in  Table  II - 2  are  theoretical; 
several  factors  influencing  the  amplitudes  of  P  waves  from  seismic  eventr 
are  not  considered  in  the  perceptibility  calculations.  These  factors  include 

•  Energy  radiation  patterns  from  the  source 

•  Anomalous  energy  attenuation 

•  Station  environment 

•  Regional  geology  in  the  vicinity  of  the  station 

Amplitudes  of  P  phases  from  the  same  event,  then,  may  vary 
considerably  at  the  same  distances,  depending  on  these  factors. 

D.  ENERGY  PROPAGATION  PATTERNS 

1.  Objectives 

The  tendency  of  magnitudes  competed  at  certain  stations  for 
events  in  the  same  areas  to<bc  consistently  higher  or  lower  than  the  average 
event  magnitude  was  observed,  and  this  suggested  the  possibility  of  determining 
distinct  energy  propagation  patterns  from  events  in  given  areas  through  in¬ 
vestigation  of  station  magnitude  residuals. 

This  investigation  was  performed  to 

•  Determine  whether  such  patterns  exist 

•  Categorize  the  patterns  found 

•  Evaluate  and  determine  the  applicability 
of  the  method 

2.  Methods  of  Investigation 

Magnitude  values  (m^)  computed  by  Texas  Instruments  (1964b) 
and  published  in  USC&GS  Earthquake  Data  Reports  (1963,  1964),  VELA 
UNIFORM  Array  Station  Bulletins  (1964,  1965),  and  LRSM  Bulletins  ( 1964,  1965) 
for  the  selected  events  are  combined  to  determir*  average  event  magnitudes. 
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AH  magnitudes  are  rounded  to  the  nearest  tenth  magnitude  unit  and,  where 
duplicate  values  are  found,  precedence  is  given  to  bulletin  or  Texas  Instru¬ 
ments  values. 


Defined  are  10  geographical  sectors  in  which  magnitude  data 
are  grouped  according  *o  record:ng- station  locations.  The  azimuthal  and 
distance  boundaries  are  different  for  Aleutian  Islands  events  and  Kurile 
Islands  events  but  correspond  to  the  same  groupings  of  stations.  The  sectors 
are  defined  as  follows: 


(1)  Canadian  Arctic  and  Greenland 

(2)  Alaska  and  Northwest  Canada 

(3)  Eastern  United  States 

(4)  Western  United  States 

(5)  Hawaiian  Islands 

(6)  Southwest  Pacific  and  Eastern  Australia 

(7)  Southeast  Asia 

(8)  India  to  Iran 

(9)  Turkey  and  Southeastern  Europe 

(10)  Western  Europe 


Magnitude  residuals  are  computed  by  subtracting  the  event 
mean  from  the  individual  magnitude  values,  and  an  average  magnitude  re¬ 
sidual  is  computed  for  each  sector  if  two  or  more  residuals  are  determined 
for  the  sector: 

n. 


P 

J 


(21) 


TT 
&  * 
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where 


_  .  ,  .  .  .th 

F.  =  mean  residual  in  j  sector 
J 

th 

t.  =  magnitude  residual  at  i  station 


n.  =  number  of  magnitudes  in  j1*1  sector 


Values  of  P.  ar'.  plotted  in  sector  boundaries  shown  in  a 
J 

polar  display  and  events  grouped  from  visual  inspection.  Values  of  P.  for 
events  associated  with  a  given  pattern  are  usually  within  ±0.  1  magnitude 
unit  of  the  mean  for  the  pattern. 

3.  Results 

Table  II- 3  lists  the  events  selected  for  the  study.  The  first 
nine  events  in  the  table  are  located  in  the  Rat  Islands  area  between  50.  3°N 
and  51.8°N  and  between  178.  1  °E  and  176.  “W.  The  Aleutian  arc  in  this 
area  trends  nearly  east-west,  with  little  apparent  curvature. 

Events  4478,  4947,  5001,  and  6536  have  similar  patterns, 
as  shown  in  Figure  II- 10.  The  pattern  is  characterized  by  relative  highs 
in  sectors  3  and  10  and  relative  lows  in  sectors  4  and  7.  The  difference  in 
mean  magnitude  residuals  between  the  Eastern  U.  S.  and  the  Western  U  S. 
averages  almost  0.6  magnitude  units  for  the  four  events  fitting  this  pattern 
( -  -  A). 

Events  2945  and  2958  are  grouped  together  as  fitting  patterr 
1-B.  This  pattern  is  characterized  by  mean  magnitude  residuals  which  are 
slightly  larger  than  0  in  the  Arctic  and  Eastern  United  States  and  slightly 
less  than  0  in  the  Western  United  States  (Figure  II- 11).  There  is  some 
evidence  that  a  nodal  plane  passes  just  off  the  West  Coast  of  the  United 
States  for  both  events  fitting  this  pattern. 
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LIST  OF  EVENTS  STUDIED 


11-40 


sclsncs  ssrvloss  division 


Figure  11-10.  Magnitude  Residual  Pattern  1-A  and  Associated  Events 
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Three  Rat  Islands  events  have  not  been  associated  with  either 
patterns  1-A  or  1-B.  Of  these  events,  2212  and  3012  show  considerable 
similarity.  In  both  cases,  there  appears  to  be  *  nodal  plane  in  a  northerly 
direction  and  a  maximum  energy  propagation  through  the  Eastern  United 
States.  The  remaining  event  (2955)  shows  some  similarity  with  pattern  1-B 
except  for  negative  mean  residuals  in  sectors  1  and  3.  Figure  11-12  shows 
the  magnitude  residual  patterns  for  events  2212,  2955,  and  3012. 

Figure  11-13  shows  the  locations  of  the  Rat  Islands  events 
studied.  The  events  fitting  pattern  1-A  are  grouped  rather  closely  on  the 
eastern  side  of  the  area  included  in  the  study.  Probably,  all  of  the  events 
are  shallow.  Events  5001  and  6536  have  depths  of  30  km  and  32  km,  respec¬ 
tively,  while  events  4478  and  4947  are  assigned  depths  of  33  km  which  implies 
a  negative  depth  encountered  in  the  hypocenter  determination. 

Events  2945  and  2958  associated  with  pattern  1-B  eu  also 
located  quite  near  each  other;  focal  depths  of  60  km  and  64  km,  respectively, 
are  determined  for  them.  Event  2955,  which  is  quite  similar  to  pattern 
1-B,  also  is  located  near  events  2945  and  2958  and  is  determined  to  have 
slightly  shallower  depth  of  50  km. 

In  this  area,  then,  similarity  in  magnitude  residual  patterns 
may  be  correlated  with  location,  dep.. .,  and  time  of  occurrence.  From 
Table  II-3,  it  can  be  observed  that  all  events  fitting  pattern  1-A  occurred 
from  August  1963  to  December  1963,  while  events  associated  with  or 
similar  to  pattern  1-B  occurred  in  March  and  April  1963. 

From  Table  II- 3  and  Figure  11-14,  it  is  evident  that  the  Kurile 
Islands  events  occurred  over  a  larger  area  than  did  the  Rat  Islands  events. 
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Rat  Islands  Events  Investi 
Classified  by  Magnitude  R 


Figure  11-13.  Location  of  Rat  Islands  Events 


Figure  11-14.  Location  of  Kurile  Islands  Events 


Four  patterns  of  magnitude  residuals  are  observed  for  the 
20  events  studied.  Pattern  19-A  shows  no  clear  indication  of  the  location  of 
a  nodal  plane,  but  all  events  fitting  this  pattern  do  apnear  to  have  a  maximum 
of  P-wave  energy  passing  through  the  Eastern  United  States.  Average  magni¬ 
tude  residuals  in  the  Western  United  States  and  Western  Europe  are  about 
equal  for  events  classif  ed  with  pattern  19-A.  Events  2098,  2319,  5422,  and 
5771  are  grouped  together  as  fitting  pattern  19-A,  as  shown  in  Figure  II - 15. 

Events  3429,  5298,  5425,  5520,  5524,  and  5578  are  classified 
as  pattern  19-B  events.  Again,  there  is  no  clear  definition  of  a  nodal  plane 
for  this  pattern,  although  residuals  at  the  westernmost  stations  in  sector  4 
indicate  the  nodal  plane  may  be  closer  to  the  U.S.  West  Coast  than  in 
pattern  19-A.  Mean  magnitude  residuals  in  sector  3  are  also  lower  than  in 
pattern  19-A,  while  the  mean  residual  in  Western  Europe  (sector  19)  is  con¬ 
siderably  higher,  indicating  that  a  maximum  of  P-wave  energy  may  be 
radiated  in  this  direction.  Figure  II- 16  shows  pattern  19-B  and  the  events 
associated  with  it. 

Events  5375,  539o,  5484,  and  5550  are  classified  as  pattern 
19-C  events.  Figure  11-17  shows  the  pattern  and  the  e’  -nts  associated  with 
it.  For  this  pattern,  there  is  considerable  evidence  of  a  nodal  plane 
passing  through  the  Western  United  States.  Also,  a  maximum  of  P-wave 
energy  appears  to  be  propagated  in  the  direction  of  the  Hawaiian  Islands. 

Mean  magnitude  residuals  in  the  Eastern  United  States  and  Western  Europe 
are  about  the  same.  It  is  interesting  to  note  that  the  difference  between 
the  mean  magnitude  residuals  in  the  Eastern  U.S.  and  the  Western  U.S. 
(sectors  3  and  4)  remains  nearly  constant  at  approximately  0.3  magnitude 
units  for  all  three  patterns.  Mean  magnitude  residuals  in  both  sectors  are 
most  positive  for  pattern  19-A  and  become  progressively  more  negative 
for  patterns  19-B  and  19-C. 
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Figure  11-15.  Magnitude  Residual  Pattern  19-A  and  Associated  Events 


Figure  11-16.  Magnitude  Residual  Pattern  19-B 
and  Associated  Events 
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Figure  11-17.  Magnitude  Residual  Pattern  19-C 
and  Associated  Events 


Events  5327  and  5523  have  markedly  similar  magnitude 

residual  patterns  and  are  grouped  together  as  pattern  19-D  (Figure  11-18). 

As  in  patteiu  19-C,  a  nodal  plane  appears  to  pass  through  or  near  the 

Aestern  United  States  (sector  4).  However;  the  maximum  of  P-wave  energy 

appears  to  be  in  the  directions  of  the  Eastern  United  States,  Western  Europe, 

and  possibly  Hawaii.  The  difference  between  mean  residuals  in  the  Eastern 

» 

U.S.  and  the  Western  U.S.  (sectors  3  and  4)  is  double  those  observed  for 
patterns  19-A,  -B,  and  -C.  It  may  be  that  a  focal  mechanism  was  operating 
for  events  5327  and  5523  that  was  different  from  that  operating  for  the  14 
events  associated  with  the  other  three  patterns. 

The  remaining  four  events  (3366,  5401,  5423,  and  5583)  fit 
none  of  the  four  patterns  observed  and  have  patterns  considerably  dissimilar 
with  one  anothei.  Observed  magnitude  residual  patterns  for  these  events 
are  shown  in  Figure  11-19. 

Little  correlation  with  event  location  and  magnitude  residual 
pattern  is  evident  for  the  Kurile  Islands  events.  Only  those  events  classified 
with  pattern  19-B  are  much  restricted  in  areal  extent.  Even  in  this  case, 
as  may  be  seen  m  Figure  11-14,  events  with  different  patterns  are  intermixed 
with  those  of  pattern  19-B  (with  respect  to  location).  Correlation  with  depth 
is  also  poor.  Depths  range  from  26  km  to  40  km  for  events  classified  with 
pattern  19-A,  from  25  km  to  80  km  for  pattern  19-B,  from  25  km  to  60  km  for 
pattern  19-C,  and  60  km  and  33  km  for  pattern  19-D  (probably  restrained 
after  a  negative  depth  was  encountered). 

Time  of  occurrence  does  not  appear  a  factor  in  separation  into 
patterns  either,  since  most  events  studied  occurred  in  October  1963.  Thus, 
while  Rat  Islands  patterns  are  separated  with  respect  to  depth  range,  location, 
and  time  of  occurrence,  no  such  correlation  can  be  made  for  Kurile  Islands 
patterns. 
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4.  Critique 

This  study  indicates  that  classification  of  events  according 
to  patterns  of  magnitude  residuals  may  be  possible.  Additionally,  t  j  study 
shows  that  there  may  be  several  such  patterns  observed  for  events  located 
in  close  proximity  to  each  other.  Also,  these  patterns  may  be  qualitatively 
correlated  with  various  focal  parameter  orientations  and  indicate  possible 
differences  in  focal  mechanisms.  However,  it  must  be  stressed  that  the  re¬ 
sults  obtained  apply  only  to  the  two  regions  studied  and  attempts  should  not 
be  made  to  extrapolate  these  results  to  other  regions. 

.  As  a  check  of  the  effectiveness  of  the  method,  magnitude 

residuals  from  the  LONGSHOT  nuclear  event  were  subjected  to  an  analysis 
similar  to  that  for  the  Rat  Islands  events.  The  observed  pattern  (Figure  11-20) 
would  not  be  classifed  as  either  1-A  or  1-B.  To  properly  assess  the  effective¬ 
ness  of  the  method,  more  events,  particularly  Aleutian  Islands  events,  need 
analysis  of  the  type  performed  iu  this  study.  Data  quality  is  probably  as  good 
as  can  be  currently  obtained.  However,  several  factors  influence  the  data: 

•  Possible  difference  in  amplitude  measurement 
criteria 

•  Errors  in  depth-distance  factors  used  in  magni¬ 
tude  computation 

•  Effects  of  station  environment  on  recorded  ampli¬ 
tudes,  affecting  computed  magnitudes 

Minimization  of  the  effe:ts  of  the  factors  just  mentioned, 
plus  considerably  more  widespread  reporting  of  short-period  P-phase 
amplitudes,  might  considerably  enhance  the  effectiveness  of  the  method 
used  in  this  study  as  a  means  of  assessing  seismic-event  focal  mechanisms 
and  might  add  a  further  criterion  for  discrimination  between  earthquakes 
and  explosions. 


11-54 


soisnos  ssrviess  division 


J 


SECTION  III 
CRITIQUE 

Data  for  this  study  were  ac  cumulated  from  nearly  every 
availab.e  source.  Of  the  highly  capable  seismograph  stawons  operating  in 
19^4,  there  were  few,  if  any,  from  which  data  were  not  collected  in  one  form 
or  another  for  this  project.  Data  collected  by  the  U.S.  Coast  and  Geodetic 
Survey,  the  International  Seismological  Centre  at  Edinburgh,  Scotland,  and 
the  Seis'iic  Data  Laboratory  at  A1 rxandria ,  Virginia,  were  made  available 
to  Texas  Instruments  Incorporated  for  this  study.  In  addition,  microfilm  re¬ 
cords  from  39  of  the  most  capable  Worldwide  Standard  Stations,  from  five 
Canadian  stations,  and  from  Matsushiro,  Japan,  were  obtained  and  subjected 
to  detailed  analyses.  These  data  were  supplemented  by  data  obtained  frorr. 
bulletins  of  various  seismograph  stations  and  networks  of  stations. 

Considerable  reliance  has  been  placed  on  statistical  proce¬ 
dures  in  the  development  of  methods  and  the  alysis  of  data.  Such  proce¬ 
dures,  ic  is  felt,  are  appropriate  for  studies  uc‘  as  this  one  which  involve 
lai-gt  quantities  of  da  a.  Also,  methods  developed  for  the  processing  of 
a- i  ay-  station  data  have,  in  this  study,  been  applied  to  single- station  data 
en. ploying  the  "world-. \r  ray"  concept.  The  two  investigations  of  such  appli¬ 
cation  have  demonstrated  the  feasibility  of  the  approach  when  applied  to 
records  originally  recorded  on  magnetic  tape. 

Conclusions  reached  in  this  study  are  based  on  well-documented 
results  usually  obtained  from  large  quantities  of  data.  Also,  these  conclusions 
and  the  results  upon  which  they  are  based  are,  in  every  case,  relevant  to  the 
original  objectives  of  the  study. 
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A  lack  of  funds  to  complete  the  processing  of  all  1964  data 
precluded  the  achievement  of  all  of  the  original  objectives;  however,  the  fol¬ 
lowing  objectives  were  satisfied: 

•  A  computer  program  was  written,  tested  and 
used  for  revision  of  January  1964  hypocenters 

•  The  program  included  several  innovations 
which  represented  improvements  in  accuracy 
and  data -handling  capability 

•  Studies  of  depth  phases  resulted  in  assess¬ 
ments  of  their  reliability  as  depth  indicators 
and  resulted  in  means  of  increasing  the  capa¬ 
bility  of  identifying  such  phases 

•  Relationships  between  magnitude  scales  were 
investigated  thoroughly,  and  tne  combination 
of  various  magnitudes  as  an  average  was  con¬ 
cluded  to  be  statistically  invalid,  even  though 
some  physical  meaning  might  be  attached  to 
such  averages 

•  P- wave  magnitud  _  computed  from  data 
within  1000  km  of  the  source  were  found  to 
be  highly  unreliable;  the  observed  variations 
were  likely  related  to  crustal  differences 

•  Patterns  were  observed  in  magnitude  residuals, 
suggesting  that  such  procedures  be  used  for 
studying  source  mechanisms  and  possibly  as 

a  discrimination  criterion 
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•  Seismic  activity  near  a  station  was  found  to 
vary  considerably  from  place  to  place 

In  addition  to  achieving  the  specific  objectives  just  outlined, 
the  study  obtained  results  which 

•  Demonstrated  the  applicability  of  array 
processing  techniques  to  single- station  data 
combined  as  a  "world  array" 

•  Established  a  means  of  easily  assessing 
epicenter-location  accuracy  from  the  dis¬ 
tribution  of  time  residuals 

•  Allowed  the  computation  of  nonlinear  or 
"exact"  confidence  regions 

Thus,  it  is  felt  that  the  original  objectives  were  met  as  well 
as  possible  within  the  financial  limitations  imposed  and,  in  some  cases, 
were  exceeded. 
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APPENDIX  A 


REVISED  HYPOCENTERS,  JANUARY  1964 


A  total  of  344  hypocenter  locations  published  by  the  USC&GS 
for  January  1964  were  input  to  the  Texas  Instruments  Incorporated  hypocenter 
program.  Of  these,  333  revisions  are  presented.  Six  locations  cannot  be 
determined  from  the  data  supplied,  and  five  revisions  yield  results  which  are 
incompatible  with  observations  of  earthquake  occurrences.  It  had  been  planned 
to  reprocess  these  ii  events,  but  circumstances  did  not  allow  it. 


The  computer  program  used  by  Texas  Instruments  (described 
in  Technical  Report  No.  Ill,  1966)  was  written  for  the  CDC  6600. 


The  following  are  the  column  headings  for  the  data  presented: 


EVENT 


SOURCE 


HR,  MIN,  SEC 

LAT. 

LONG. 

DEPTH 


N 


Chronological  number  assigned  to  each 
event  located  by  the  USC&GS 

TI  and  PDE  for  Texas  Instruments 
revised  hypocenter  and  input  prelim¬ 
inary  hypocenters  located  by  the  USC&GS, 
respectively 

Origin  time  (GCT)  in  hours,  minutes, 
seconds,  and  tenths  of  seconds 

Geographical  latitude  in  degrees 

Geographical  longitude  in  degrees 

Depth  of  earthquake  foci  in  km;  (*depth  re 
strained  to  33  km  after  negative  depth  en¬ 
countered  iteration  process) 

Number  of  stations  used  for  each  location 


SD 

A,  B,  and  C 


MB  and  MS 

PEG 


Standard  deviation  for  each  location  in 
seconds  (shown  for  TI  data  only) 

Semiminor  axis,  semimajor  axis,  and 
inclination  of  the  7  0-percent  confidence 
ellipses  (shown  for  TI  data  only) 

Magnitudes  computed  on  m^  and  M  scales, 
respectively 

Seismic  region  based  on  Gutenberg  and 
Richter's  (1954)  region  as  modified  by 
Flinn  and  Engdahl  (1964) 
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APPENDIX  B 


ANALYSIS  OF  HYPOCENTER  REVISION 

The  most  striking  single  result  of  the  hypocenter  revision 
performed  by  Texas  Instruments  Incorporated  is  the  difference  in  depth 
distribution  of  che  located  hypocenters.  Texas  Instruments  results  indicate 
that  less  than  60  percent  of  the  earthquakes  located  in  January  1964  occurred 
at  depths  of  70  km  and  shallower  and  that  less  than  50  percent  occurred  at 
depths  of  50  km  and  shallower.  The  preliminary  results  published  by  USC&GS, 
however,  indicate  that  more  than  two-thirds  of  the  earthquakes  occurred  at 
depths  of  70  km  and  shallower  and  more  than  60  percent  occurred  at  depths 
of  50  km  and  snallower.  This  result  is  primarily  due  to  the  decrease  in  the 
number  of  events  with  depths  restrained  to  33  km  by  Texas  Instruments. 

As  shown  in  Table  B-l,  151  of  the  USC&GS  hypocenters  were 
restrained  to  depths  of  33  km  (about  45  percent  of  the  total)  while  TI  results 
show  only  43  events  so  restrained.  Of  the  108  events  restrained  to  33  km 
by  USC&GS  but  not  by  Texas  Instruments,  33  depths  were  placed  shallower 
than  33  km  and  7  5  deeper  than  33  km.  These  changes  in  depths  resulted  in 
47  events  placed  at  depths  greater  than  50  km  and  35  events  at  depths  greater 
than  70  km  rather  than  at  33  km. 

Since  the  USC&GS  usually  places  the  accuracy  of  its  published 
depth  determination  at  ±25  km,  this  figure  is  used  to  establish  the  number  of 
depths  determined  by  Texas  Instruments,  which  differs  significantly  from 
that  determined  by  the  USC&.GS.  Using  this  criterion,  109  of  the  TI  depth 
determinations  are  found  to  differ  appreciably  from  the  USC&GS  depths. 

The  results  discussed  here  should  not  be  construed  to  imply 
that  109  of  the  USC&GS  depth  determinations  are  incorrect.  Due  to  lack  of 
time  and  funds,  Texas  Instruments  results  could  not  be  subjected  to  as  rigorous 
and  extensive  quality  control  as  the  USC&GS  routinely  applies  to  their  data. 
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DEPTH  DISTRIBUTION 


Depth 

(km) 

N  umber 

TI 

o i  Events 

CGS 

1  to  10 

17 

1 1  to  20 

14 

12 

21  to  30 

18 

5 

31  to  40 

•A. 

T 

*j! 

76 

166 

41  to  50 

30 

27 

51  to  60 

20 

10 

61  to  70 

14 

9 

71  to  100 

35 

20 

101  to  150 

42 

25 

151  to  200 

26 

17 

201  to  250 

11 

12 

251  to  300 

1 

1 

301  to  350 

3 

3 

351  to  400 

4 

1 

401  to  450 

4 

6 

451  to  500 

3 

2 

501  to  550 

5 

7 

551  to  600 

6 

8 

601  to  650 

_ 2_ 

2 

Total 

333 

333 

43  events  restrained  to  33-km  depth 
151  events  restrained  to  33-km  depth 
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To  say  that  one  set  of  results  is  'better  1  than  another  is  a  rather  dubious 
statement.  To  further  illustrate  this,  a  comparison  with  the  International 
Seismological  Center  (ISC)  results  shows  both  agreement  and  disparity  with 
USCScGS  and  Texas  Instruments  results.  The  observed  differences  in  depth 
determination  illustrate  the  extent  of  the  depth -determination  problem.  Con¬ 
siderable  study  of  this  problem  is  strongly  recommended.  Data  of  all  types 
should  be  used,  and  overconfidence  in  computer  processing  and  precon¬ 
ceived  ideas  should  be  avoided  in  such  studies. 

The  significance  of  epicenter  shifts  is  somewhat  more  difficult 
to  assess.  Only  90  of  the  333  epicenters  are  unchanged  in  the  revision  process. 
A  total  of  109  epicenters  are  shifted  0.  1°  latitude  and/or  longitude  and  134  epi¬ 
center  shifts  exceed  0.  lc  latitude  and/or  longitude. 

Another  indication  of  the  degree  of  changes  in  epicenter  loca¬ 
tions  is  that  165  of  the  TI  epicenters  do  not  include  the  USC&GS  location 
within  the  70-percent  confidence  region.  In  other  words,  probably  about 
115  of  the  TI  epicenters  do  differ  from  USC&GS  epicenters. 

When  changes  are  considered  in  either  epicenter  location  and 
focal  depth,  or  both,  appreciabiechanges  are  found  to  have  been  made  in  half 
(a  total  of  168)  of  the  input  hypocenters  processed.  Criteria  for  significance 
are  epicenter-location  changes  exceeding  0.  1°  latitude  and/or  longitude  and 
focal -depth  changes  exceeding  2  5  km,  or  both. 

In  evaluating  the  differences  in  the  results,  two  factors  must 
be  considered:  differences  in  the  methods  and  the  data  used.  With  regard 
to  method,  the  Texas  Instruments  hypocenter  program  differs  in  two  princi¬ 
pal  ways  from  the  USC&GS  program: 

•  Travel-time  corrections  are  applied  to  observed 
times  for  nearly  100  stations,  and  these  corrections 
depend  on  areas  in  which  the  events  occurred 

•  Larger  time  residuals  are  allowed  than  in  the 
USC&GS  program 
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Table  B-2 

AVERAGE  NUMBER  OF  STATIONS  IN  EACH  QUADRANT  FOR  VARIOUS  SEISMIC  REGIONS 


The  first  of  these  two  differences  probably  has  the  greater  effect  on  hypo- 
cencer  determinations. 

Data  used  by  Texas  Instruments  differ  in  both  quantity  and 
quality.  In  locating  333  iiypocente rs,  TI  used  more  than  9000  station  readings 
compared  to  4270  used  by  the  USC&GS.  Of  course,  the  use  of  moie  data 
does  not  necessarily  improve  the  accuracy  of  hypocenter  determinations.  An 
important  factor  is  the  distribution  of  these  data.  Table  B-2  compares  quad¬ 
rant  distribution  of  data  used  by  Texas  Instruments  to  determine  hypocenter 
locations  in  various  seismic  regions.  Quadrant  1  refers  to  epicen^er-to-station 
azimuths  ranging  from  0°  to  90  J;  quadrant  2,  from  90°  to  1  8 0 c ;  quadrant 
3,  from  180°  to  270°;  and  quadrant  4,  from  270°  to  360°. 

Based  on  the  average  number  of  stations  used  in  each  quadrant, 
the  quadrant  distribution  of  data  used  by  Texas  Instruments  is  somewhat  im¬ 
proved  over  that  used  by  the  USC&GS  in  regions  1,  5,  7,  8,  10,  15,  16,  17, 

18,  19,  29,  30,  34,  and  44.  The  proportions  of  additional  data  used  for  hypo¬ 
center  determinations  in  regions  20  through  26  and  in  48  which  are  in  the 
fourth,  or  northwest,  nuadrant  —  result  in  considerably  greater  weight  on  data 
to  the  northwest  in  the  TL  determinations  than  in  the  USC&GS  determinations. 

Con. paring  the  quantity  of  data  on  a  regional  basis,  it  may  be 
seen  that  TI  uses  approximately  three  times  as  much  data  as  does  the  USC&GS 
for  determinations  in  regions  17,  18,  19,  20,  21,  22,  25,  26,  29,  and  30; 
approximately  twice  as  much  data  is  used  for  determinations  in  regions  1,  7, 

22,  and  48.  At  least  40  percent  more  data  is  used  in  the  regions  listed  in 
Table  B-2,  except  in  regions  34  and  44.  Thus,  the  quantity  and/or  quadrant 
distributions  of  data  are  considerably  enhanced  in  nearly  all  regions. 

The  quality,  or  a^  racy,  of  arrival  times  of  most  data  used 
was  probably  somewhat  improved  over  that  supplied  to  the  TJ3C&GS  since,  in 
a  considerable  number  of  cases,  data  analysis  was  performed  after  the  USC&GS 
preliminary  hypocenter  locations  were  published. 
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In  summary  .appreciable  changes  in  about  half  of  the  preliminary 
hypocenturs  for  January  1964  are  effected  through  use  of  the  Texas  Instruments 
hypocenter  program,  A  portion  of  this  result  stems  from  the  apnlication  of 
station/source-region  travel-time  corrections  and  another  portion  from  the 
use  of  considerably  more  data,  which  are  generally  better  distributed  and 
somewhat  improved  in  overall  quality  o~  accuracy.  Probably  the  most  signif¬ 
icant  results  are  the  sharp  decrease  in  the  number  of  events  with  depths  re¬ 
strained  to  ^3  km  and  the  increase  in  the  proportion  of  events  with  interme¬ 
diate  focal  depths.  However,  these  results  require  further  investigation 
as  does  the  depth  problem  in  general. 

Further  improvement  in  hypocenter  determination  accuracy 

requires 

•  Considerably  more  knowledge  of  variations  in 
travel  times  resulting  from  effects  near  both 
source  and  station 

•  Availability  of  more  high-quality  data  (particu¬ 
larly  in  the  Southern  Hemisphere) 

•  Development  of  methods  for  weighting  data  on 
the  basis  of  quality  and  distribution 

It  is  also  recommended  that  alternatives  to  the  least-mean- 
squares  method  of  hypocenter  determination  be  explored. 
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APPENDIX  C 


LOCAL  SEISMICITY 

The  number  and  density  of  seismograph  stations  in  the  Middle 
East  and  in  Southern  Asia  precludes  the  location  of  many  earthquakes  of  small 
or  even  moderate  magnitudes.  Consequently,  in  the  area  south  of  the  USSR 
and  China,  most  seismicity  evaluations  are  based  on  relatively  large  events. 

In  this  report,  an  alternative  method  has  been  used.  Rather 
than  evaluating  the  seismicity  of  the  entire  region,  the  seismicity  of  an  area 
included  by  a  radius  of  1000  km  around  the  more  capable  seismograph  stations 
has  been  evaluated.  The  following  are  the  procedures: 

•  Determine  the  perceptibility  limits  for  values  of  m^ 
from  2.  0  to  5.  5  ir  increments  of  0.  5  units 

•  Determine  the  epicentral  distances  and  magnitudes 
of  all  events  within  1000  km  of  the  station 

•  Count  the  number  of  events  with  magnitudes  of  2.  0 
±  0.  25,  2.  5  ±  0.  25,  etc.  ,  which  occur  within  the 
calculated  limits  of  perceptibility  for  m^  =  2.0,  2.5,  etc. 

•  Normalize  this  number  to  a  radius  of  1000  km  in  each 

6  2  2 

case;  i.e.  ,  multiply  by  10  km  and  divide  by  r  where 
r  is  the  calculated  perceptib;lity  limit 

•  Plot  the  logarithms  (base  10)  of  these  numbers  vs 
magnitude 

•  Fit  a  straight  line  to  the  plotted  points  by  the  least- 
squares 

Data  obtained  from  the  analyses  of  short-period  vertical  seis¬ 
mograms  recorded  in  1963  at  Istanbul,  Turkey;  Shiraz,  Iran;  Quetta,  Pakistan; 
and  Shillong,  India  were  used  to  accomplish  the  study.  Results  are  shown  in 
Figures  C-l  through  C-4. 
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NUMBER  OF  EARTHQUAKES 


science 


On  the  basis  of  the  results  obtained,  the  rates  of  earthquake 
recurrence  in  the  vicinities  of  Istanbul  and  Quetta  are  nearly  the  same  and 
only  slightly  larger  near  Shillong.  However,  the  recurrence  rate  for  earth¬ 
quakes  near  Shiraz  is  nearly  twice  that  for  Quetta.  Thus,  in  1963,  the  num¬ 
bers  of  earthquakes  with  magnitudes  of  m^  =  5.  0  ±  0.  25  within  1000  km  of 
the  stations  are 


Istanbul 

Shiraz 

Quetta 

Shillong 

At  magnitudes  of  m  = 

b 

numbers  are 

Istanbul 

Shiraz 

Quetta 

Shillong 


13 

1 

8 

3 

.  0  ±  0.25,  however,  the  estimated 

263 

309 

148 

108 


Thus,  if  seismicity  were  based  on  only  numbers  of  events, 
the  Shiraz  area  would  be  considered  the  most  seismic;  however,  since  more 
larger  events  occur  near  Istanbul,  the  total  seismic  energy  release  in  this 
area  is  highest.  The  conclusion  reached  is  that  the  area  surrounding  Istanbul 
is  the  most  seismic  of  those  studied. 
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APPENDIX  D 


WORLDWIDE  ARRAY  PROCESSING 
A.  INTRODUCTION  AND  SUMMARY  OF  RESULTS 

Digital  processing  of  records  from  stations  forming  a  world¬ 
wide  array  (as  reported  by  Texas  Instruments  Incorporated  in  Special  Report 
No.  4)  was  initially  undertaken  as  a  possible  means  of  enhancing  depth-phase 
identification.  The  present  study  is  a  continuation  of  the  investigation  into 
the  problem  of  energy  propagation  and  magnitude  det_rmination.  In  the  pre¬ 
vious  study,  errors  may  have  been  present  in  the  data  because  of  hand-digitizing 
from  papei  records;  this  difficulty  has  been  overcome  in  the  present  study 
since  data  were  recorded  on  magnetic  tape  and  computer-digitized. 

Worldwide  array  processing  techniques  have  been  applied  to 
data  from  a  selected  suite  of  five  Kurile  Islands  events  and  the  LONGSHOT 
nuclear  explosion  in  the  Aleutian  Islands.  Table  D-l  presents  the  U.S.  Coast 
and  Geodetic  Survey  'n  nocenter  information  on  these  six  events. 


Table  D-l 

EVENTS  ANALYZED  BY  WORLDWIDE  ARRAY  PROCESSING 


Date 

H 

Lat. 

Long. 

h 

mb 

No. 

Sta. 

Location 

Assigned 

No. 

1/15/64 

02  23  47.4 

45. 3N 

150.  6E 

45 

5.3 

10 

Kurile  Islands 

15 

3/16/66 

08  44  32.  8 

41.  8N 

146. 8E 

140 

5.  7 

12 

Kurile  Islands 

76 

7/1/64 

09  46  49.6 

45.  2N 

150.  3E 

75 

4.  8 

10 

Kurile  T.s lands 

183 

11/12/64 

05  16  26.  2 

47.  IN 

146.  6E 

328 

4.  8 

8 

Kurile  Islands 

317 

12/17/64 

05  18  34. 8 

45.  4N 

150. IE 

17 

5.  3 

10 

Kurile  Islands 

352 

10/26/65 

21  00  00.  1 

51. 4N 

179.  2E 

0 

6.0 

33 

Aleutian  Islands 

LS 
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The  processing  sequence  shown  in  Figure  D-l  was  applied  to 
all  data.  Then,  the  processed  data  were  studied  for  information  relevant  to 
four  areas  of  special  interest. 

1.  Depth-Phase  Identification 

Short-period  vertical  tiaces  from  each  station  were  gathered 
and  composited  to  form  a  single  record  with  individual,  equalized  traces 
displayed  side  by  side  in  the  order  of  increasing  epicentral  distance;  41  pP 
phases  could  be  identified  on  the  composite  record  as  compared  to  12  pP 
phases  reported  by  the  individual  stations.  Stacked  traces  from  three  of  the 
five  earthquake  events  clearly  showed  a  pP  phase. 

2.  Digital  Filtering  Techniques 

In  previous  worldwide  array  processing,  mos:  filtering  was 
analog  filtering  accomplished  electrically  on  playback.  Digital  filtering  is 
generally  superior  to  analog  filtering  because  digital  filters  are  stronger  and 
more  selective  and  the  filtered  output  can  be  input  into  further  processing. 
Development  of  digital  filtering  capabilities,  therefore,  was  considered  im¬ 
portant  to  further  progress  in  worldwide  array  processing. 

Each  of  10  filters  designed  was  applied  to  16  selected  tract  , 
and  the  results  were  displayed  side  by  side  on  a  single  recoid.  Results  of 
the  study  were  twofold: 

•  A  library  of  10  assorted  filters  is  stored 
on  magnetic  tape  for  future  use 

•  There  is  visual  evidence  of  each  filter's 
effectiveness  on  various  types  of  traces 
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3.  Power  Density  Spectrum  Studies 

The  power  density  spectrum  of  the  short-period  vertical  trace 
was  computed  for  each  station.  The  differences  in  power  spectra  recorded 
by  different  stations  for  the  same  event  appear  to  depend  only  on  the  local 
crustal  structure  at  the  stations.  Power  spectrum  results,  therefore,  have 

% 

been  used  to  compute  magnitude  corrections  for  nine  stations. 

4.  Energy  Attenuation 

Energy  attenuation  curves  were  plotted  from  average  P  energy 
in  the  first  minute  after  the  P-phase  arrival  at  10  stations  recording  LONG- 
SHOT  at  nearly  the  same  azimuth.  Comparison  of  these  results  with  pres¬ 
ently  used  Q  curves  shows  discrepancies  as  large  as  one  m^  unit;  as  a  result, 
the  need  for  Q-curve  revision  is  obvious.  The  results  of  this  study  also  in¬ 
dicate  that  energy  attenuation,  when  average  power  is  considered,  may  not 
be  as  irregular  as  would  be  expected  considering  the  large  standard  deviations 
obtained  for  magnitudes  calculated  from  a  single  cycle. 

-2Kr 

0 

The  energy-vs-distance  curve  may  be  approximated  by  — 

r 

for  A  >  56°.  For  A  <  56°,  the  increasing  energy  with  increasing  distance  be¬ 
havior  needs  further  investigation;  the  cause  may  possibly  be  diffraction 
effects  at  the  Mohorovicic  discontinuity. 

B.  DIGITAL  PROCESSING  TECHNIQUES 

All  data  used  in  this  study  were  originally  recorded  in  analog 

form  on  magnetic  tape  by  Long-Range  Seismic  Measurement  Vans.  Digitizing 

* 

was  done  by  a  Texas  Instruments  DARC  (Data  Analysis  and  Reduction  Com¬ 
puter)  directly  from  composite  analog  magnetic  tapes  supplied  by  the  Seismic 
nata  Laboratory  in  Alexandria,  Virginia.  This  procedure  eliminated  any 


❖ 
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errors  which  might  have  been  present  4  »  previous  studies  (Special  Report 
No.  4)  as  a  result  of  hand-digitizing.  The  digitized  records  were  in  a  format 

)Je 

compatible  with  TIA.C  (Texas  Instruments  Automatic  Computer),  which  was 
used  to  accomplish  all  data  processing. 

A  TIAC  trace-gather  ’-outine  transferred  the  short-period 
vertical  trace  from  each  station  onto  a  composite  record  for  each  event. 

Other  routines  were  then  used  to  edit  the  composite  records.  Editing  included 
the  shifting  of  traces  to  line  up  the  first  peak  of  the  P-arrival  trace,  zeroing 
to  remove  spikes  or  other  bad  segments  from  the  traces,  and  trace  polarity 
reversal  to  make  the  direction  of  the  first  motion  of  all  of  the  traces  the  same. 

The  records  were  then  equalized  over  the  first  minute  after 
the  P  arrival  by 

•  Computing  the  average  trace  amplitude  for 
the  first  minute 

•  Obtaining  the  ratio  of  the  desired  average 
amplitude  to  the  computed  average  amplitude 

•  Multiplying  every  point  on  the  original  trace 
by  this  amplitude  ratio  and  output 

The  output  records  from  equalization  have  the  seme  average 
amplitude  and,  consequently,  the  same  power  over  the  first  minute. 

Power  spectra  were  computed  by  taking  the  Fourier  transform 

of  the  autocorrelation  of  the  trace: 

T 

max 

autocorrelation  =  ^  ^  tr  (t)  •  tr  (t  -T)  (1) 

t  =  T 

o 
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where 


A  T 

tr{t) 

tr(t-r) 


0,  At,  2At,  ....  t 

max 

sample  rate  of  the  correlation 
value  of  trace  at  time  t 
value  of  trace  at  t-r 


Therefore,  power  spectrum 


uu 

-f 


autocorrelation 


“if  T  j 
e  dr 


(2) 


The  term  "stacking"  refers  to  a  simple  algebraic  sum  of  N 

traces.  The  theory  behind  stacking  is  that  the  noise  in  each  of  the  N  traces 

is  different.  Therefore,  the  total  noi9e  in  all  traces  is  randomly  out  of  phase 

and  tends  to  cancel  in  the  algebraic  sum.  Signals  such  as  pP  which  arrive  at 

all  stations  at  a  nearly  constant  interval  from  the  P  phase,  however,  tend  to 

be  in  phase  on  all  traces  and,  therefore,  will  add  in  the  sum.  The  pP  ampli- 

1/2 

tude-to-noise  ratio  on  the  stacked  trace,  then,  is  improved  by  a  factor  of  N 

The  digital  filters  designed  in  the  study  were  computed  by  TIAC. 
The  computer  reads  in  the  parameters  describing  the  desired  filter  and  out¬ 
puts  the  filter  onto  magnetic  tape  ready  for  use.  This  filter  is  fed  into  the 
computer  against  the  data  trace,  resulting  in  a  filtered  output  trace.  The 
main  advantage  of  a  digital  filter  over  an  analog  filter  is  that  the  output  from 
the  digital  filtering  process  can  be  input  into  later  processing  steps.  Another 
advantage  is  that  the  digital  filter  can  be  designed  to  be  much  "sharper"  than 
the  analog  filter.  Figure  D-2  shows  typical  power  spectra  for  both  types  of 
filters. 

C.  DEPTH-PHASE  IDENTIFICATION 

The  application  of  worldwide  array  processing  techniques  to 
depth-phase  identification  has  been  reported  in  Special  Report  No.  4.  The 
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Figure  D-2.  Power  Spectra  of  Typical  Analog  and  Digital  Filte 


present  effort  has  limited  the  work  in  this  area  to  two  processes  to  avoid 
duplication;  i.  e.  , 

•  Composite  records  are  edited,  equalized, 
and  displayed  in  the  order  of  increasing 
epicentral  distance  on  a  single  record 

•  All  of  the  traces  for  each  event  aligned  on 
the  P  arrival  are  stacked  together,  form¬ 
ing  a  single  trace 

Figures  D-3  through  D-8  show  results  of  the  first  process. 
Depth-phase  identification  is  considerably  improved  by  forming  these  arrays. 
Table  D-2  lists  the  number  of  traces  on  which  pP  and  PcP  can  be  identified 
from  the  array  as  compared  to  the  number  of  these  phases  reported  by  in¬ 
dividual  stations. 


Table  D-2 

DEPTH-PHASE  IDENTIFICATION- 
SINGLE  TRACE  VS  WORLDWIDE  ARRAY 


E 

V 

e 

n 

t 

PP 

PcP 

Individual 

Trace 

Array 

Individual 

T  race 

Array 

15 

1 

9 

0 

0 

76 

6 

7 

0 

0 

183 

3 

9 

0 

0 

317 

1 

6 

0 

1 

352 

1 

10 

0 

1 

LS 

0 

0 

28 

21 
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Figure  D-3.  Earthquake  Event  15 
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Figure  D-4.  Earthquake  Event  76 
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Figure  D-6.  Earthquake  Event  317 
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Figure  D-7.  Earthquake  Event  352 
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Figure  D-8.  Array  Ensemble  Presentation  of  LONGSHOT  Recording  (2  of  2) 
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The  second  process  results  in  positive  pP  identification  for 
three  of  the  five  processed  earthquake  events  (15,  76,  and  352),  which  indi¬ 
cates  t:\at  stacking  may  be  developed  into  a  useful  tool  for  depth-phase  iden¬ 
tification.  Figure  D-9  shows  the  stacked  traces  for  the  three  events  of  interest. 

One  interesting  discrepancy  appeared  after  study  of  the  stacked 
traces  of  event  352.  The  stacked  trace  indicates  a  pP-P  interval  of  14.  75 
sec  which,  at  an  average  distance  of  50°,  indicates  a  depth  of  approximately 
60  km  as  compared  with  17  km  reported  by  the  USC&GS.  For  events  15  and 
76,  the  stacked  traces  indicate  depths  of  45  km  and  140  km,  respectively, 
which  exactly  agree  with  the  USC&tGS  value. 

D.  DIGITAL  FILTERING  STUDIES 

Since  previous  array  processing  did  not  utilize  digital  Liters 
to  any  extent,  as  assortment  of  high-pass,  low-pass,  and  bandpass  digital 
filters  with  frequency  ranges  as  listed  in  Table  D-3  were  designed  and  applied 
to  selected  traces  for  the  purpose  of  better  determining  their  usefulness  in 
data  processing.  Selected  from  the  data  were  16  representative  traces  — 
two  from  each  earthquake  event  and  six  from  LONGSHOT.  All  10  filters 
were  applied  to  each  of  the  selected  traces.  For  each  input  trace,  a  single 
output  record  was  received,  with  the  10  filtered  outputs  of  the  same  input 
trace  displayed  side  by  side. 

These  filtered  outputs  show  considerable  improvement  over 
the  raw  data.  Here,  filtering  is  definitely  superior  to  the  previous  analog 
filtering.  Figure  D-10  shows  two  examples. 

As  a  result  of  this  study,  a  library  of  10  filters  is  stored  on 
magnetic  tape  ready  for  immediate  use  in  the  future.  In  addition,  there  is 
evidence  of  the  effectiveness  of  each  of  the  filters  on  various  types  of  traces. 
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EVENT  15  pP-P  INTERVAL  =  12.25  SEC 


P  pP 


EVENT  352  pP-P  INTERVAL  s  14. 75  SEC 


P 


Figure  D-9.  Stacked  Traces  Showing  Depth  Phases 


Table  D-3 
DIGITAL  FILTERS 


Filter  Type 

T 

(sec) 

high-pass 

0-1.5 

high-pass 

0  -  1.0 

high-pass 

0  -  2.  0 

bandpas  s 

0.  5  -  1.  5 

bandpass 

1.  0  -  1.  5 

bandpass 

0.  67-2.0 

bandpass 

1.0  -  2.0 

bandpas  s 

0.  5  -  1. 0 

low -pass 

0.  5  -  00 

low  -  pas  s 

1.  0  -  00 
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FILTERED  TRACE  (.67<t<2.0) 


The  linear  position  of  10  stations  along  a  single  azimuth  from 
the  LONGSHOT  event  (Figure  D-ll)  provides  an  opportunity  to  study  energy 
attenuation  with  distance.  The  short-period  vertical  traces  from  each  of  the 
10  stations  selected  for  this  study  are  equalized  over  a  1-min  time  window, 
beginning  with  the  first  P  arrival. 

The  vertical  component  of  the  total  energy  arriving  at  the 
station  in  the  first  minute  is 
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VIE  =  total  energy  in  the  first  minute  of  the 
equalized  trace 


G 

IR 


=  instrument  gain 

=  amplitude  of  calibration  signal  (peak  to  peak) 
equivalent  ground  motion  given  in  daily  tape  logs 


=  instrument  response  of  f 
Section  II 


pmax 


as  calculated  in 


EQS  =  equalization  scalar 

=  amplitude  of  chosen  peak  on  equalized  trace 
amplitude  of  same  peak  on  unequalized  trace 


As  a  result  of  equalization,  E  is  a  constant  and  can  be 

VI 

omitted  from  Equation  (3).  The  vertical  component  of  the  energy,  then,  is 
given  by 

E y  a  [g2  x  (IR)2  x  (EQS)2]  (4) 


Figure  D-12  shows  results  of  the  calculation  of  the  vertical 
component  of  P  energy  for  the  10  selected  stations  as  a  function  of  distance. 
Assuming  *hat  energy  radiates  spherically  from  the  source,  the  energy  per 
unit  area  perpendicular  to  the  direction  of  wave  propagation  would  diminish 
as  follows: 


E(r ) 


E 

source 


x 


-2K 


e 

2 


r 


r 


(5) 


The  vertical  component  of  the  energy  at  distance  r  is 


•  Z'.'. 

E  (r)  =  cos2  (e  )  E(r)  =  E  cos2  (e  )  e  r  (6) 

v  o  source _ o _ 

2 

r 

where  e^  is  the  angle  that  the  direction  of  propagation  makes  with  the  normal 

to  the  earth's  surface.  The  energy-vs-distance  curve  in  Figure  D-12  closely 

-2Kr 

follows  the  expected  —  curve  for  £  >  56°. 

r  2 
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Figure  D-ll.  Recording  Stations  tor  LONGSHOT  Event 
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For  small  distances,  the  cos  (e  )  term  would  he  expected  to 

c 

decrease  vortical  energy  with  decreasing  distances;  however,  the  effect 
would  be  slight  and  could  not  account  for  the  large  decrease  observed. 


Further  study  into  this  matter  is  needed;  one  possible  cause 
might  be  diffraction  at  the  Mohorovicic  discontinuity  (Figure  D-13).  This 
approach  also  explains  the  anomaly  in  the  Rocky  Mountain  area. 

Figure  D-14  shows  the  spherical  divergence  effect  and  the 
possible  diffraction  effect  which,  in  combination,  resemble  tha  LONGSHOT 
results.  The  fact  that  these  results  are  in  the  form  of  the  simple  functions 
predicted  by  making  the  assumption  of  spherical  divergence,  which  is  equiva¬ 
lent  to  the  assumption  of  homogeneous  crust  and  mantle,  indicates  that  one 
could  expect  more  consistent  energy  or  magnitude  calculations  by  considering 
the  average  P-wave  amplitude  rather  than  the  amplitude  of  the  maximum  cycle. 

Figure  D-12  has  been  transformed  into  a  more  familiar  repre¬ 
sentation  called  the  Q  curve.  The  unified  magnitude  is  given  by 

mb  =  log  T  +  ° 


where  Q  is  an  empirical  correction  factor.  The  Q  curve  is  the  p’ot  of  this 

correction  factor  as  a  function  of  distance  and  depth.  The  commonly  used 

Q  curve  for  surface-focus  events  is  compared  with  the  Q  curve  calculated 

from  LONGSHOT  data  (Figure  D-15).  The  units  on  the  vertical  axis  and  the 

relative  positions  of  the  two  curves  are  arbitrary.  The  data  indicate  that 

present  Q  curves  may  introduce  errors  as  large  as  one  magnitude  unit  into 

the  calculation  of  m,  .  It  seems,  therefore,  that  a  Q-curve  revision  is  needed 

b 

for  improvement  of  the  scale. 

The  presence  of  anomalies  in  the  Q  curves,  such  as  the  one 
seen  in  the  LONGSHOT  curve,  are  very  dependent  on  source  location,  there¬ 
fore,  it  may  be  necessary  to  compute  sets  of  Q  curves  for  every  source  region. 
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Figure  D-14.  Two  Effects  Limiting  Energy  Arriving  at  Surface  of  Earth 
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Figure  D-15.  Q  Curves  (Surface  Focus) 
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An  important  result  of  this  study  is  that  a  technique  has  been 
developed  to  apply  data  processing  techniques  to  the  precise  computation  of 
Q  curves,  but  the  method  is  now  somewhat  limited  because  of  lack  of  data 
recorded  on  magnetic  tape. 

F.  POWER  DENSITY  SPECTRA  STUDIES 

Power  density  spectra  are  computed  for  all  stations  and  all 
events.  The  spectra  are  computed  over  a  1-min  time  window,  beginning  with 
the  P  arrival. 

Figure  D-16  presents  power  density  spectra  for  all  stations 
recording  event  183,  and  Figure  D-17  shows  power  density  spectra  for 
LONGSHOT  recorded  at  stations  lying  along  a  single  azimuth.  For  all  six 
events,  the  power  density  spectra  seem  to  have  the  same  general  appearance. 
Deviations  from  this  general  appearance  do  not  seem  to  vary  in  a  systematic 
fashion  as  distance  increases.  LeBlancand  Howell  (1966)  also  observed  this 
result. 

The  frequency  of  maximum  power  is  calculated  from  the  power 
density  spectrum  of  each  station  event  and  is  subtracted  from  the  frequency 
of  maximum  power  for  each  station  recording  each  event.  Figure  D-18  pre¬ 
sents  the  results  of  this  calculation  where  the  frequency  of  maximum  power 
shows  no  well-behaved  distance  dependence. 

Figures  D-19  and  D-20  give  magnitude  and  depth  dependence 
of  the  average  frequency  of  maximum  power  for  the  six  events  studied.  Event 
focal  depth  seems  to  have  no  effect  on  energy  frequency,  although  Moham- 
madioum  (1965a  and  1965b)  has  observed  this  effect  for  Alaskan  earthquakes. 
A  tendency  toward  higher  frequencies  for  larger  events  is  indicated. 
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Figure  D- 16.  Power  Density  Spectra  for  All  Stations  Recording  Event  183 
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Figure  D-18.  All  Events 
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Figure  D-19.  f  Vs  Magnitude 
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Figure  D-20.  f  Vs  Focal  Depth 
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Results  of  TI  studies  indicate  that  the  power  density  spectrum 
of  the  energy,  as  it  travels  through  the  earth,  remains  essentially  unchanged 
over  the  distances  studied.  The  differences  in  power  spectra  between  dif¬ 
ferent  stations  for  a  given  event  would  arise,  therefore,  from,  the  filtering 
effect  of  the  complex  geology  of  the  crust  at  the  station.  If  this  is  true,  the 
power  spectra  from  several  events  occ  rring  in  the  same  region  and  re¬ 
corded  at  a  given  station  should  be  similar.  Figure  D-21  presents  the  simi¬ 
larity  in  power  spectra  for  three  KurUe  Islands  events  recorded  at  TFO  and 
BMO. 

For  every  station,  then,  a  crustal  operator  G  may  be  defined; 
it  operates  on  the  pov.or  spectrum  $  of  the  energy  arriving  at  the  bottom  of 
the  crust  to  give  the  power  spectrum  recorded  at  station  i.  In  general, 

$  =  Mf)  $  i  Ma)  (3) 

and 

*.  (f,  e  ,  0)  =  G.(f,  e  ,  0)  Hi)  (9) 

i  o  xo 

If  all  events  considered  are  from  the  same  region, 

G.  (f,  eo,  0)  a  G.  (f)  s  g.(f)  (10) 

Equation  (9)  for  the  case  where  events  are  in  the  same  region 

reduces  to 

*.'(f)  =  g-(f)  *(f)  (H) 

Another  subscript  is  added  to  Equation  (11)  to  identify  tht 

event: 


$.'.(f)  =  g.  (f)  $  (f) 


(12) 
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Equation  (1Z)  represents  the  crustal  operator  for  station  i 
operating  on  the  power  density  spectrum  of  event  j,  yielding  the  recorded 
power  density  spectrum  at  station  i  for  event  j. 

The  power  density  spectrum  data  may  he  used  from  n  events 
to  compute  the  average  crustal  operator  for  station  i,  according  to  the  fol¬ 
lowing  equation: 


i  £ 


$:.(f) 

JU _ 

$.<f) 

j 


(13) 


A  third  subscript  added  to  Equation  (12)  to  identify  source 

region  gives 


Vjk<«  ■  %<«  vf)  (14) 

Solving  for  the  average  crustal  operator  for  station  i,  given 
power  density  spectrum  data  from  n  events, 


g 


$.(f) 

J 


1 

as  — 
n 


n 


L 

j=i 


$ 


1  m  ,  •> 

~  y),. 

m  L,  ij 
i=l 


(15) 


The  average  crustal  operator  g.,  (f),  as  given  in  Equation  (15), 

liC 

can  be  used  to  compute  the  portion  of  the  energy  arriving  at  the  base  of  the 
crust  at  station  i  from  source  region  k,  which  will  be  filtered  cut. 

The  average  crustal  operator,  then,  is  a  direct  indicator  of  the 
errors  introduced  into  magnitude  calcu'adons  based  on  the  energy  arriving 
at  station  i  and  can  be  used  to  calculate  a  station  magnitude  correction  factor. 
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Power  density  spectrum  data  from  the  six  events  processed  in 
this  study  are  used  to  calculate  g.^(f)  for  nine  stations.  The  results  are  pre¬ 
sented  in  Figure  D-22.  The  magnitude  calculated  from  a  cycle  of  frequency 
f  recorded  at  one  of  these  nine  stations  may  be  corrected  for  station  effects 
by  being  changed  by  an  amount  corresponding  to  the  value  of  the  magnitude 

correction  at  f  for  that  station. 

c 

G.  RECOMMENDATIONS 

The  following  recommendations  for  further  study  are  based  on 
results  of  the  present  study: 

•  Digital  processing  techniques  are  shown  to 
be  of  value  in  many  areas  of  interest  to  seis¬ 
mology.  Further  investigation  of  these  areas 
is  severely  limited  because  of  lack  of  data  re¬ 
corded  on  magnetic  tape,  so  a  worldwide  en¬ 
semble  of  stations  equipped  with  magnetic  tape 
recording  instruments  is  recoi  'mended;  such 
instrumentation  would  make  dig  fal  processing 
techniques  readily  available  for  seismological 
research 

•  Based  on  discrepancies  between  LONGS*iOT 
energy  attenuation  and  the  commonly  accepted 

Q  curves,  revision  of  Q  curves  is  recommended 

•  The  average  energy  in  the  P  packet,  at  least 
for  LONGSHOT,  was  shown  to  be  a  reasonably 
well-behaved  predictable  function  of  distance; 
therefore,  it  is  the  expressed  opinion  that  a 
magnitude  scale  based  on  average  P  energy 
would  be  superior  to  present  scales  based  on 
the  amplitude  of  a  single  cycle 
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REPORTS 

A.  SEMIANNUAL  TECHNICAL  REPORTS 

1.  No.  I,  28  April  1965 

The  purpose  of  this  first  report  is  to  introduce  and  establish 
ue  background  for  discussions  to  be  presented  in  future  reports.  Since  studies 
have  shown  there  are  errors  in  determining  epicenters  and  depths  of  focus  of 
seismic  events,  this  report  discusses  development  ox  techniques  to  be  studied 
for  improvement  in  determining  epicenters,  depths  of  focus,  magnitudes,  and 
annual  number  of  se’smic  events,  However,  hypocenter  determination  is 
emphasized.  A  discussion  of  the  data  necessary  for  the  study  is  also  included. 

2.  No.  II,  28  October  1965 

This  is  a  progress  report  concerning  the  determination  of 
seismic  event  parameters  such  as  size,  number  of  occurrences,  and  accurate 
locations.  Specifically,  however,  this  report  presents  the  accomplishments 
under  Contract  C- 104-65  following  the  first  semiannual  report  and  observa¬ 
tions  concerning  preliminary  findings.  Since  the  first  report  emphasizes 
determination,  the  bulk  of  this  second  report  deals  with  hypocenter  deter¬ 
mination.  Discussions  of  magnitude  studies,  seismicity  evaluation,  and  two 
special  studies  of  depth  phases  are  also  included. 

3.  No.  JTI,  28  April  1966 

This  technical  report  shows  the  progress  achieved  toward 
accomplishing  the  objectives  initiated  in  Contract  C- 104-6  5,  reports  the 
status  of  various  studies,  and  st..:es  plans  for  completing  the  contract. 

Section  II  describes  the  data  used  and  the  status  of  data  receipt  and  data  re¬ 
duction  and  handling.  Section  III  treats  the  hypocenter  determination  program 
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and  associated  studies.  Section  IV  discusses  1963  magnitude  data  investiga¬ 
tion  results  and  their  impact  on  magnitude  determinations  for  1964  events,  as 
well  as  the  computer  program  to  be  used  in  calculating  magnitudes  for  1964 
events.  Section  V  discusses  future  plans  and  scheoales  for  completing  the 
current  contract.  The  various  topics  detailed  in  the  bod/  of  this  report  are 
summarized  in  Section  I. 

Be  SPECIAL  REPORTS 

1.  No.  1,  Revaluation  of  Seismicity  for  I960  and  1963 

Included  are  revisions  and  corrections  to  previously  published 
data.  Data  are  prepared  for  comparison  with  1964  seismicity  as  well  as  for 
combination  with  1964  data  to  permit  further  statistical  analysis.  Interpre¬ 
tations  of  the  1960  and  1963  data  are  modified,  revised,  and  supplemented. 

2.  No.  2,  Correlation  of  Time  Residuals  with  Magnitude 

This  report  presents  the  results  of  an  investigation  to  determine 
whether  an  empirical  correlation  may  be  observed  between  time  residuals 
and  magnitudes  on  a  station/ sou  rce- region  basis.  Results  of  a  study  con¬ 
cerning  the  value  of  the  average  time  residual  as  a  function  of  increasing 
magnitude  are  given.  Included  is  a  study  to  determine  whether  station/source- 
region  time  correlations  expressed  as  a  function  of  magnitude  will  increase 
the  accuracy  of  hypocenter  determinations. 

3.  No.  3,  Depth  Phases 

The  objectives  of  this  report  are  to  appraise  the  likelihood  of 
identifying  depth  phases  from  paper  seismograms  and  to  evaluate  differences 
in  the  reliability  of  various  depth- indicator  phases. 

4.  No.  4,  Results  of  Worldwide  Array  Feasibility  Investigation 

'  The  object  of  this  report  is  to  enhrr.ee  the  identification  of 
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depth  phases  by  combining  sirgle-trace  recordings  of  short-period  vertical 
seismograms  (USC&GS  Worldwide  Standard  Stations)  in  ensemble  displays  and 
applying  some  multichannel  computer  processing  techniques  to  the  ensembles. 

5.  No.  5,  Method  for  Treating  Cumulative  Errors  in  Epicenter  Deter¬ 
minations 

The  method  which  treats  cumulative  errors  and  time- residual 
patterns  is  studied  to  determine  a  theoretical  epicenter  shift  toward  the  true 
location.  The  method  is  applied  to  data  from  LONGSHOT  —  a  surface  event 
with  a  computed  depth  of  57  km. 
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